{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import math\n",
    "import matplotlib.pyplot as plt\n",
    "from scipy.fftpack import fft,ifft\n",
    "import pandas as pd\n",
    "from scipy.optimize import minimize"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoMAAAHnCAYAAADD41dfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3XuUZWdZJ+DfSwICgi2YgE6Ck7jCcIeAEVEHbVAwjCQgogOCYBACCijOoIQlI+IddTkKAsolhICCSeQWCMYR6DA4yJCEJgQDs0JUaKIrQaS5LMIlvPNHnSZFparrnO5zqXP286xVK7X32Wfvt9LJqV+/3/72V90dAACG6SaLLgAAgMURBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAE7ctEFLJOjjjqqjzvuuEWXAQCwrUsuueST3X30dscJgxM47rjjcvHFFy+6DACAbVXVP49znGFiAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGGS57N699gUATIUwCAAwYEcuugA4qI1dwL17k/37b7x/z545FQQAq0VnkOWza9eiKwCAlaEzyM62seN3oCOoEwgAU6EzCAAwYMIgAMCAGSZmuRgeBoCp0hlk8Tw7EAAWRhgEABgwYRAAYMDcM8j8jfsg6cQ9ggAwYzqD7AweJA0AC6EzyPx5kDQA7Bg6gwAAA6YzyOLpCALAwugMMhueHQgAS0EYBAAYMMPEHJ6tun8eFwMAS0FnkNnxuBgA2PF0Bjk8W3X5PC4GAJaCziAAwIDpDDIbOoIAsBR0BgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYZHu7d9+w1jAAsFKEQQCAAbM2MTfYqvu3d2+yf//mr1uDGACWms4g49m1a9EVAAAzoDPIDbbq8h3oCOoCAsDK0RkEABgwnUG2pyMIACtLZxAAYMCEQQCAARMGWU0elA0AYxEGAQAGzAQSVsPGLuBWD8o2GQYAvo7OIKvLg7IBYFs6g6yGjR0/D8oGgLHoDAIADJgwCAAwYIaJWU2GhwFgLDqDAAADJgwCAAyYMMhwWaUEAIRBAIAhEwYBAAbMbGKGZf2w8GZL1pmFDMDA6AwybJasA2DgdAYZlvWdP0vWAYDOIADAkOkMMlw6ggCgMwgAMGTCIADAgAmDcDBWKQFgxQmDAAADJgzC4dA5BGDJmU0M620MdlYpAWDF6QzCdg5nlRKdQwB2OJ1BWG9j12+zVUqsbwzACtEZhMNlfWMAlpjOIExqu/WND9Y51DUEYIcRBuFgphHedA4B2MGEQZi27TqHALCDCINwOIQ8AJacCSQAAAMmDAIADJhhYpglw8gA7HA6gwAAAyYMDoml0QCADYRBAIABc8/gqjlY52+zdXTXc3/bfHkGIQA7gM7g0FgNAwBYR2dw1Rysy6QTBQBsoDMIADBgOoMwT+vv19x4D6eOLQALIAwOibCx87iHE4AFEwZhntYHch1BAHYA9wwCAAyYMAgAMGCGiWFRDA8DsAPoDAIADJgwCAAwYMIgAMCACYOwE+3e/fUPqAaAGREGAQAGTBiEZaJjCMCUebQM7BTWLQZgAXQGYaeybjEAc6AzCDvFwdYtPrCtYwjAlOkMwrLRMQRginQGYRls7BDqCAIwJTqDAAADpjMIO5HOHwBzojMIADBgOoOwTHQMAZgynUEAgAETBmGVWb4OgG0Ig7AKhD4ADpEwCAAwYMIgAMCAmU0Mq8haxgCMSRiEZbZV6FvPWsYAHIQwCKtifeizljEAYxIGYZkJfQAcJhNIAAAGTBgEABgww8SwygwbA7ANYRBWwbRCn3sPAQbHMDEAwIAJgwAAAyYMwhDt3r35A6oBGBz3DAKbr2TivkGAQdAZBL6e5esABkVnELCSCcCAHVIYrKpvTHJdd18/5XqAWdpsODgR/gAGbKxh4qq6SVX9ZFW9taquSfLhJP9SVR+qqt+vqjvOtkxg6gwHA5DxO4PvTPK3SZ6d5PLu/mqSVNVtkzwgye9W1Ru6+zWzKROYCsPBAGwwbhj8oe7+8sad3f2pJH+V5K+q6qZTrQyYP+EQYHDGGiZeHwSr6lnbHQOsEM8kBFhp23YGq+qc9ZtJTkzy/JlVBADA3IwzTPyZ7n7igY2qeskM6wHmwXAwACPjDBP/1obtX5lFIQAAzN+2YbC7/zFJquqo0fanZl0UsMO5jxBgZUzy0Okzk5w6q0KAHcZ6xQCDMMnaxDWzKoCdzQOqAVbWJJ3BnlkVjMeDgpknD6gGGIRJwqDOIAydtY0BVs4kYfDZM6uCGxzspvyNv4A345cy82LoGGAljBUGq+rO3X35rIuZp6p6eJIfSXK7JC/q7r9ZcEnj8QuYRTJ0DLByxu0MXlBVe5L8Wnd/bIb1jKWqzkzy0CTXdPfd1+0/OckfJzkiycu7+3e3Okd3vzHJG6vqNkn+IMnOCIMH++XqFzAAMGXjhsE7J3lykouq6k1Jfqu7r51dWds6K8mfJDn7wI6qOiLJi5I8KMm+JO+rqjdnLRj+zob3P6G7rxl9/5zR+4DN+MsHwEob69Ey3f2l7n5hkrtkLWi9t6p+vapuPdPqtq7nXUk2Pvz6vkmu7O6ruvtLSV6X5GHd/cHufuiGr2tqzfOTvK27L93qWlV1elVdXFUXX3vtIvMvAMD0TfKcwXT3dd39B0nukeS6JJdW1TNnUtnkjkny8XXb+0b7tvL0JD+U5JFV9ZStDurul3b3Sd190tFHHz2dSg/Vnj26NOxsViYBWDqTzCZOVR2XtSHjOyX59iSfTfLbWbvnbtE2e/TNls9G7O4XJHnB7MqBFeYvJQArY9zZxJclOTbJx5J8OMkVSd6RtXvtPjKz6iazL8kd1m0fm+TqBdUCALAUxu0M/miSq7p7J69C8r4kd6yq45N8IsmjkvzkYksCANjZxp1A8tHu7qq6Y1W9oqr+ZNaFHUxVvTbJe5Lcqar2VdXPdPdXkjwtyYVZ61ye090fWmSdAAA73UT3DCZ5dZLnJXl+klTV3ZP8cnc/btqFHUx3P3qL/RckuWCetQCxTB3AEptoNnGSm3T325JcnySjVUnufvC3AINilRyApTJpZ/Dq0T15nSRVVUluMfWqgOVimTqApTVpGHxGkpcn+daqOi3JyUlWas1iAIAhmSgMdvc/jdb/fXiSeyW5KMmZsygMAIDZm/Sh07/Z3c9Jct7o68CawAAALKFJJ5AcU1Vfe3ZfVR2d5G+nWxKwUixRB7CjTXrP4JOTXFhVV2ZtEskrkzxr6lUBy2mSiSMmmwDsCOMuR3d2kkuTvD/JU5P8RZKvJHl4d185u/IAAJilcYeJXzU69glZC4LHJfn3JI+tqkfOpjQAAGZtrM5gd789ydsPbFfVkUnumrUZxffLaDIJAADLZdJ7BpMko3WALxt9vXqqFQGrY/duS9QB7HCHFAYBJrJ+ibrN1jEWEAEWRhgEZmfPnht3BNc/ZsY6xgALJwwC82UdY4AdZdKHTqeqPFcQAGBFbNsZrKpz1m8mOTHJ82dWEQAAczPOMPFnuvuJBzaq6iUzrAcAgDkaJwz+1obtX5lFIcCKskQdwI627T2D3f2PSVJVR422PzXrogAAmI9JZhOfmeTUWRUCDIzuH8COMMls4ppZFQAALMQkYbBnVgUAAAsxyTCxziAwG5aoA1iYSTqDz55ZFQAHbLdE3e7dX7+kHQCHZezOYHdfXlUXJTmluz9TVU9JcvMkL+7uL82swh2gqk5JcsoJJ5yw6FJgNVmiDmBhJl2O7ptHQfA7kzwpyW2SvGz6Ze0s3X1+d5++a7uOBQDAkpnknsEk+XJVHZnkcUme393nVNXFM6gLAIA5mDQMviDJB7I2PHzGaN+tploRAABzM1EY7O6zq+r1Sa7v7i9U1QlJ3jOb0gAAmLVJ7xlMd3+uu78w+v7K7j5t+mUBbGL37rXHz1x0kRnFAFMy6TAxwOyMO4vYZC6AqREGgeWxZ4/HzwBM2cTDxAAArI5tw2BVPaiqXlZVJ462T599WQAAzMM4w8Q/l+S0JM+pqtsmOXG2JQEAMC/jDBNf292f7u5nJnlwku+acU0AAMzJOGHwrQe+6e4zkpw9u3IA1tm92yNkAGZs22Hi7n5TklTVUd39ye5+4ezLAtiCWcQAUzXJbOIzZ1YFAAALMUkYrJlVAQDAQkwSBntmVQAAsBA6gwAAAzbJcnTPnlkVABvt3p3s3Zvs328JOoAZGrsz2N2XJ0lV3XR25QBssGvXoisAWGmTdAZTVS9P8oiq+nySq5NcluQyj5sBpm7PHh1BgDmYKAwmuX+S23f3l6vqmCT3SnLP6ZcFAMA8TBoG/z7JbZJc092fSPKJJBdMvSoAAOZiktnESfLSJBdV1TOr6v5V5WYeYHlY3g7gRiYNg69Jck7WOoo/l+T/VNVHp14VAABzMekw8b7ufu76HVX1DVOsB+DQmXACMLFJw+DeqvqF7v7jAzu6+4tTrglgjVAHMHOThsHbJ/mhqnpWkkuTfCDJ3u4+d+qVAUzD+nsENz7E+gChExiwicJgd/9E8rWh4bsluUeS704iDALLwUOsAb7OpA+dPinJFd39+ax1Bi+dSVUA49jY4dus87e+6+eeQoAbmXQ28dlJrj+wUVVHVdVDp1sSwGHQ+QOYyKT3DF7X3dcd2OjuT1bVryd5y3TLAhjDxg6fzh/AxCbtDF5VVQ/ZsO9m0yoGAID5mrQz+PQkb6uqn8ra0nR3S+Kh0wAAS2rS2cT/UlXfmeThSU5MclmS/zaLwnaSqjolySknnHDCoksBDofhY4Abqe5edA1L46STTuqLL7540WUAs+S+Q2BFVNUl3X3SdsdNes8gAAArRBgEABiwicJgVV1UVd80+v4pVfWMqjKbGABgSU3aGfzm7v7MaBLJk5LcJsnLpl8WwBLavfvGq6IA7HCTPlrmy1V1ZJLHJXl+d59TVWZUAMttfYDbuKSdiSTAips0DL4gyQeS3DzJGaN9t5pqRQCLZkk7YEAmfc7g2VX1+iTXd/cXquqEJO+ZTWkAc7K++6cjCAzMpJ3BdPfn1n1/ZZLTploRAABzM1EYrKoHJnlMkk8nuTxrK5Bc3t1fnEFtAADM2KSdwdckeerofffM2rJ0d0tinTZgmEw+AZbcpGHwyu5+w+j7c6ddDMDSM/kEWDKThsGLquoXk/xRW9QYWEWTdvNMPgGW3KRh8G5J7p7kWVV1SZK9SfZ2ty4hAMASmvTRMo9Ikqq6RW4IhveLIWMAgKU06Wzi2yb5xSS3S/IPSc7u7rNmUBcAAHMw6drEr0vy2STnJ7llkndX1X2nXhXATmTtYWAFTXrP4Ld19++Nvn9LVf1lkr/I2lAxAABLZtIw+Kmqumd3X5Yk3X1VVd1yBnUBLB+ziIElNGkYfHKS86rqfyf5YJK7Jvno1KsCAGAuJr1n8AtJ7pPknVmbRPKBJI+edlEAAMzHpJ3BN3T3fZKcc2BHVd0vyd9PtSqAneJgy80lhoaBpTdWZ7CqfqKqfjfJravqLlV1xLqXXzqb0gB2oENdbs5MZGCHGrcz+HdJbp7kiUn+MMmdqurTSa7O2tAxwGqa13JzlrIDFmSsMNjdn0hydlV9tLv/LvnaA6iPT/LhGdYHAMAMTXrP4BVV9Rv5+hVIPj/9sgAAmAcrkAAADJgVSAAABswKJACzst1jadbb6nUTSoAZO9wVSO4WK5AAQ3G4wWy7x9Ic6mNrAA7DRGGwuz9cVfdJ8vCsLUW3N8l/n0VhAEtvksfSeLQMsCAThcGqemCSxyT5dJLLk1yWpGdQFwAAczDpMPFrkjx19L57Zq1DeLckJ0y5LgAA5mDSMHhld79h9P250y4GAID5Gndt4rOr6hlJ3lNV7hEEAFgR4z50+lWjY781yU9V1T9X1Zur6jeq6sdnVx4AALM07trEb0/y9gPbVXVk1mYT3yvJfWPIGODgtpslbBYxsCBjhcGqelySP8xad/D8JE8bPXj6shnWBgDAjI07TPyrSR6U5M5JPpbkt2dWEQAAczNuGPxMd7+/u6/p7v+RtaHhwaiqU6rqpfv37190KQAAUzVuGPy2qjq9qu5fVUcnueksi9ppuvv87j59l6WiAIAVM+5zBp+btYdMPybJPZLcqqouSPKBJJd192tnVB8AADM07mzil67frqpjsxYO75HkvyQRBgEAltC4s4mru7+2BnF370uyL8kFWx0DAMDON+49g++sqqdX1bev31lVN6uqB1bVq5I8fvrlAQAwS+PeM3hykickeW1VHZ/k00lukbUw+TdJ/md3751NiQAAzMq49wxel+TFSV5cVTdNclSSL3T3p2dZHAAAszVuZ/BruvvLSf5lBrUAADBn294zWFUPqqqXVdWJo+3TZ18WAJvavXvta1bHA4MzTmfw55KcluQ5VXXbJCfOtiQAAOZlnNnE13b3p7v7mUkenOS7ZlwTALBousqDMU4YfOuBb7r7jCRnz64cAIBNCKczs20Y7O43VdWd122/cLYlAQAwL+POJr6gqvYk+bXu/tgM6wHgUGzVMdm7N9m/f/PX9+w5+Lm2eh1YKeOGwTsneXKSi6rqTUl+q7uvnV1ZACS5cYg7WLjbyq5d06wIWDHjPnT6S0leWFUvS/K0JO+tqtck+f3u/uwsCwRgg83CnS4fcIgmeuj0aCWSP6iqlyT5hSSXVtWfdfcfzKQ6gKHbGOKEO7YzrUkWh9KF3oz/Vne8cWYTf01VHVdVJyd5YpJvT/LZJL89i8IAWHJmfy4/txgMwlidwaq6LMmxST6W5MNJrkjyjiQvSvKRmVUHAExmWp04XejBGHeY+EeTXNXdPctiAJijQ5mBnAgHsGLGnUDy0SSpqjsmOSPJF7r7abMsDIAFMjzItLiHccebaAJJklcneV6S5ydJVd09yS939+OmXRgAU3CwX3xmILNs/CVlJiYNgzfp7rdV1W8nSXdfPgqEAMyDgMaycQ/jjjdpGLy6qo5P0klSVZXkFlOvCoCda9xhukmG9fyCh4WZNAw+I8nLk3xrVZ2W5OQkl0+9KgBWg2E92PEmDYNPzFoAfHiSeyW5KMmrpl0UADvYuF08w3rLzZ/bYEwaBo9J8uPd/dok51XV7ZK8JckDpl4ZAAAzN2kYfHKSC6vqo0m+muSsJM+adlEALJCOEAzKuCuQnJ3k0iTvT/LUJH+R5CtJHt7dV86uPAAAZmnctYlfNTr2CVkLgscl+fckj62qR86mNAAAZm3cFUjenuTtB7ar6sgkd83aJJL7JTlvJtUBACRuX5ihSe8ZTJJ091eSXDb6evVUKwJgNfjlDUth3GFiAABWkDAIADBgwiAAwIAJgwAAAyYMAgAMmDAIADBgwiAAwIAJgwAAAzb4MFhVd6mqP62q86rqZxddDwA7yO7da1+wwuYeBqvqm0fB68NVdUVVfc8hnufMqrqmqi7f5LWTq+ojVXVlVZ1xsPN09xXd/ZQkP5HkpEOpBQBgWS2iM/jHSf66u++ctbWNr1j/YlXdrqpuvWHfCZuc56wkJ2/cWVVHJHlRkodkbf3kR1fVXavqHlX1lg1ftxu959Qk78669ZcBYCF0I5mzuYbBqvqmJN+f5BVJ0t1f6u5PbzjsB5K8qapuPnrPk5K8YOO5uvtdST61yWXum+TK7r6qu7+U5HVJHtbdH+zuh274umZ0rjd39/cmecyUflQAgKVw5Jyv9x1Jrk3yyqq6V5JLkvxCd3/+wAHdfW5VHZ/kdVV1bpInJHnQBNc4JsnH123vS/LdWx1cVbuTPCLJNyS5YItjTklyygknbNagBABYXvMeJj4yyX2SvKS7753k80ludE9fd/9ekuuSvCTJqd39uQmuUZvs660O7u493f3z3f3k7n7RFsec392n79q1a4IyAJgZQ6kwNfPuDO5Lsq+73zvaPi+bhMGqun+Suyd5Q5LnJnnahNe4w7rtY5NcfUjVArA6DiU87t2b7N8/+Xv37Jn8WrAgc+0Mdve/Jvl4Vd1ptOsHk/zD+mOq6t5JXpbkYUlOS3LbqvrNCS7zviR3rKrjq+pmSR6V5M2HXTwAw2RUiBU3785gkjw9yZ+PgtpVWQt8690yyY9390eTpKoen+SnN56kql6bZHeSo6pqX5LndvcruvsrVfW0JBcmOSLJmd39oVn9MAAsiUPp1h3oCOr0scLmHga7e28O8jy/7v67DdtfzlqncONxjz7IOS7IFpNBAGDmDud+xkMdmj5AcGVCg1+BBAB2HEPTzNEihokB4OtN2gU7lO7ZPDtmh3MtQ9PMmc4gAMtJ9wymQmcQgMWbtAumewZTozMIADBgOoMAsBWdRwZAZxAAYMB0BgFgJ9GNZM50BgEABkxnEIDlo3sGU6MzCAAwYMIgAMCACYMAAAMmDAIADJgwCAAwYMIgAMCACYMAAAMmDAIADJgwCAAwYMIgAMCACYMAAAMmDAIADJgwCAAwYMIgAMCACYMAAAMmDAIADJgwCAAwYMIgAMCACYMAAAMmDAIADJgwCAAwYMIgAMCACYMAAAMmDAIADJgwCAAwYMIgAMCACYMAAAMmDAIADJgwCAAwYMIgAMCACYMAAAMmDAIADJgwCAAwYMIgAMCACYMAAAMmDAIADJgwCAAwYMLgGKrqlKp66f79+xddCgDAVAmDY+ju87v79F27di26FACAqRIGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABgwYRAAYBK7d699rQhhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABO3LRBQAAzM3u3Yd/jr17k/37D/1ce/Ycfg1TpDMIADCpXbsWXcHU6AwCAMMxja7cgY7gDuvwHSqdQQCAARMGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABiwIxddAADAUtmzZ9EVTJXOIADAgAmDAAADJgwCAAyYMAgAMGDCIADAgAmDAAADJgwCAAyYMAgAMGCDD4NVdZeq+tOqOq+qfnbR9QAAzNNCwmBVHVFV76+qtxzGOc6sqmuq6vJNXju5qj5SVVdW1RkHO093X9HdT0nyE0lOOtR6AACW0aI6g7+Q5IrNXqiq21XVrTfsO2GTQ89KcvIm7z8iyYuSPCTJXZM8uqruWlX3qKq3bPi63eg9pyZ5d5K3H84PBQCwbOYeBqvq2CQ/kuTlWxzyA0neVFU3Hx3/pCQv2HhQd78ryac2ef99k1zZ3Vd195eSvC7Jw7r7g9390A1f14zO9ebu/t4kjznsHxAAYIkcuYBr/lGSX05y681e7O5zq+r4JK+rqnOTPCHJgyY4/zFJPr5ue1+S797q4KraneQRSb4hyQVbHHNKklNOOGGzBiUAwPKaa2ewqh6a5JruvuRgx3X37yW5LslLkpza3Z+b5DKbnfIg19rT3T/f3U/u7hdtccz53X36rl27JigDAGDnm/cw8fclObWq/ilrw7cPrKrXbDyoqu6f5O5J3pDkuRNeY1+SO6zbPjbJ1YdULQDAiptrGOzuZ3f3sd19XJJHJXlHdz92/TFVde8kL0vysCSnJbltVf3mBJd5X5I7VtXxVXWz0XXePJUfAABgxezE5wzeMsmPd/dHu/urSR6f5J83HlRVr03yniR3qqp9VfUzSdLdX0nytCQXZm3G8jnd/aG5VQ8AsESqe8vb6djgpJNO6osvvnjRZQAAbKuqLunubZ+hvBM7gwAAzInO4ASq6tpsMmR9mHYl2T/lc87Couuc5/Vnea1pnnsa5zrccxyV5JOHWQPTtej/V+dtWX7eRde5Cp+h0z7vED5D/2N3H73dQcLgglXVS7v79EXXsZ1F1znP68/yWtM89zTOdbjnqKqLxxmCYH4W/f/qvC3Lz7voOlfhM3Ta5/UZegPDxIt3/qILGNOi65zn9Wd5rWmeexrnWvSfK9M3tD/TZfl5F13nKnyGTvu8PkNHdAaBQ7ZT/lYLsIx2ymeoziBwOF666AIAltiO+AzVGQQAGDCdQQCAARMGAQAGTBgEABgwYRCYiqr6xqp6VVW9rKoes+h6AJZJVX1HVb2iqs6b97WFQWBLVXVmVV1TVZdv2H9yVX2kqq6sqjNGux+R5LzuflKSU+deLMAOM8lnaHdf1d0/s4g6hUHgYM5KcvL6HVV1RJIXJXlIkrsmeXRV3TXJsUk+Pjrs+jnWCLBTnZXxP0MXRhgEttTd70ryqQ2775vkytHfYr+U5HVJHpZkX9YCYeKzBWDSz9CF8YENTOqY3NABTNZC4DFJXp/kx6rqJVmRJZoAZmDTz9Cq+paq+tMk966qZ8+zoCPneTFgJdQm+7q7P5/ktHkXA7BktvoM/bckT5l3MYnOIDC5fUnusG772CRXL6gWgGWz4z5DhUFgUu9LcseqOr6qbpbkUUnevOCaAJbFjvsMFQaBLVXVa5O8J8mdqmpfVf1Md38lydOSXJjkiiTndPeHFlknwE60LJ+h1d2LvD4AAAukMwgAMGDCIADAgAmDAAADJgwCAAyYMAgAMGDCIADAgAmDwOBV1fVVtXfd1xmLrumAqjqvqr5j9P3nJnjfnqr64Q37nlFVL66qo6vqr6ddK7CcrE0MkHyhu0+c5gmr6sjRw2UP5xx3S3JEd191CG9/bdZWNrhw3b5HJfml7r62qv6lqr6vu//ucGoElp/OIMAWquqfqup5VXVpVX2wqu482v+NVXVmVb2vqt5fVQ8b7f/pqjq3qs5P8jdVdZNRJ+5DVfWWqrqgqh5ZVT9YVW9Yd50HVdXrNynhMUnetEldR1XVe6rqR0bbvzSq5bIjXwcCAAACz0lEQVSqet7osPOSPLSqvmF0zHFJ/kOSd49ef+Po/MDACYMAyS02DBP/13WvfbK775PkJUmeOdr3K0ne0d3fleQBSX6/qr5x9Nr3JHl8dz8wySOSHJfkHkmeOHotSd6R5C5VdfRo+7Qkr9ykru9Lcsn6HVV1+yRvTfKr3f3WqnpwkjsmuW+SE5N8Z1V9f3f/W5L/m+Tk0VsfleQv+4Zlpy5Ocv8x//0AK8wwMcDBh4kPdOwuyVq4S5IHJzm1qg6Ew5sn+fbR9/+ruz81+v4/Jzm3u7+a5F+r6p1J0t1dVa9O8tiqemXWQuLjNrn2tyW5dt32TZO8PclTu/uidbU8OMn7R9u3ylo4fFduGCp+0+ifT1h3rmuy1ikEBk4YBDi4L47+eX1u+MysJD/W3R9Zf2BVfXeSz6/fdZDzvjLJ+Umuy1pg3Oz+wi9kLWge8JWshdIfTnIgDFaS3+nuP9vk/W9M8odVdZ8kt+juS9e9dvPR+YGBM0wMMLkLkzy9qipJqureWxz37iQ/Nrp38PZJdh94obuvTnJ1kuckOWuL91+R5IR125217t6d1814vjDJE6rqVqNajqmq242u8bkke5KcmbUu4Xr/Kcnl2/2gwOrTGQQY3TO4bvuvu/tgj5f5jSR/lOSyUSD8pyQP3eS4v0ryg1kLXf8vyXuT7F/3+p8nObq7/2GL67w1awHybw/s6O7rq+pRSc6vqs9094ur6i5J3jPKpp9L8tisDQMnayHw9VkbJl7vAaPzAwNXN9xLDMC0VdWtuvtzVfUtWZvQ8X3d/a+j1/4kyfu7+xVbvPcWSd45es/1U67rXUke1t3/Ps3zAstHGASYoarak+Sbk9wsye9191mj/Zdk7f7CB3X3Fw/y/h9OckV3f2yKNR2dtYD5xmmdE1hewiAAwICZQAIAMGDCIADAgAmDAAADJgwCAAyYMAgAMGDCIADAgP1/VB0dXFC507oAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 720x576 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "data00 = pd.read_csv(\"energy_2000.csv\")  \n",
    "x00 = data00['x']\n",
    "y00 = data00['y']\n",
    "xerr00 = data00['xerr']\n",
    "yerr00 = data00['yerr']\n",
    "\n",
    "\n",
    "plt.figure(figsize=(10,8))\n",
    "plt.loglog()\n",
    "plt.errorbar(x00, y00, yerr=yerr00, xerr=xerr00, fmt=',r')\n",
    "plt.xlabel(\"Energy (keV)\")\n",
    "plt.ylabel(r\"$keV^2 (Photons{\\ }cm^{-2}{\\ } s^{-1}{\\ } keV^{-1})$\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoMAAAHnCAYAAADD41dfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3X20rWdZH+rfDQGhopvSBPWANDhCgQCyxS3S49EsVDCeI0KxUhSrBQGpwCkdhw5gSOuotlo9Do9HRXpAAwQrNAQRAlHaIntTLbUkYRGjgZ5A/YjRAZSyEQbhI979Y66VzCzW2mvOteb3e11j7LHnfOfXvT72nL99P+/zPNXdAQBgmO6y7AIAAFgeYRAAYMCEQQCAARMGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABiw85ZdwDo5//zz+8ILL1x2GQAAh7r22ms/2t0XHHY/YXAKF154Ya655ppllwEAcKiq+uNJ7meYGABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABO2/ZBTAQW1tfeGx7Ozl7NrnkkjsfP316ERUBANEZZNlOnFh2BQAwaDqDLMZ+3b7dbqFOIAAsjTDI8giBALB0hokBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGGS2trb2X2AaAFhJlpbh6M61q8je2ywjAwArSWeQ2bOrCACsDZ1Bjk63DwDWns4gq895iAAwNzqDrJ69wc95iAAwNzqDrAfnIQLAXOgMsnp0/ABgYXQGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABgwYRAAYMCEQQCAATtv2QWwIra2vvDY9nZy9mxyySX7P+b06XlWBAAsgM4g53bixLIrAADmSGeQEV0+ABgknUEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAzlt2AczY1tb+x7e3k7Nnk0su2f/206fnVREAsMJ0BofkxIllVwAArBidwU2jw3duu51T3ycASCIMsun2DpvvDpePHxcMARgww8QMj+FyALidziCbTdcPAM5JGGR9HTQEPD5jWhgEgHMyTMxmOe4Q8NbWwcvzAMAG0hlkfe3t+h11pvB4+DPBBICBEQbZHLMKbSaYADAgwiDo/AEwYM4ZBAAYMGEQAGDAhEEAgAFzziAcZr/ZxtYyBGBD6AzCtMw2BmCD6AzCYcY7f0ddyxAAVpQwCNMQAgHYMIaJAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBiEWdrauvMi1QCw4iwtA8e13w4l48csRwPACtMZhFmzQwkAa0RnEI7LDiUArDFhEGZJCARgzRgmBgAYMGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYsPOWXQBsvK2t0d/b28nZs8kll9xx2+nTy6gIAG6nMwiLdOLEsisAgDvRGYR52+3+7XYIdQMBWCHCICyKEAjACjJMDAAwYMIgAMCAGSaGZTpoprHzDAFYEGEQVsX4TOO9IXH3eiIYAjBTwiAs00EdwPHwZzkaAOZIGFwnhgw3196fqZ8xAAsiDK6a8Y7QXvsNGe4lRAAAUzCbeN0YMgQAZkhncNXo7AEAC6QzCOtqa+vcpwwAwAR0BmHdHLTsjK4yAEegMwjrzDmkAByTziCsG7uTADBDwiCsKyEQgBkwTAwAMGDCIADAgAmDAAADJgwCAAyYMAgAMGDCIADAgAmDAAADJgwCAAyYMAgAMGDCIGyara07tqoDgEPYjg42xW4A3N5Ozp4dXbdlHQCHOFJnsKq+uKruOutigBk5cWLZFQCwJibqDFbVXZI8NcnTknxdks8k+aKq+kiSq5O8orv//7lVCRxutwu42yHUFQRgApMOE78zyX9I8pIkN3T3XyVJVd0nyWOT/KuqelN3/+p8ygQmJgQCMIVJw+C3dvfn9h7s7o8leWOSN1bV3WZaGQAAczfROYPjQbCqXnTYfQAAWA+Hdgar6orxq0lOJvmpuVUEAMDCTDJM/Inufubulap6+RzrAQBggSYZJv6Xe67/yDwKAebIQtQAHODQzmB3/7ckqarzu/ujO5NGgHWw30LUiRnHANxumkWnL5tbFcD8WYgagH1Msx1dza0KYD4sRA3AIaYJgz23KoD5EgIBOMA0w8Q6gwAAG2aaMPiSuVUBAMBSTBQGq+oh3X3DvIsBAGCxJu0MXl1Vl1XVA+ZaDbBY1h8EGLxJw+BDkrw3yZmq+rmqumCONQEAsCAThcHu/mx3/0KShya5OcnvVdWPVdWXzLU6YD52O4Lb28mZM7qDAAM2zQSSdPet3f0zSR6R5NYk11XVC+dSGbAYFqMGGLRp1hlMVV2Y0ZDxg5M8IMlfJvmJJD8z68KAObIYNQA7JgqDVXV9kvsn+ZMk709yY5LfTvKyJB+YW3XAfAmBAIM3aWfw7yT5UHfbhQQAYINMOoHkg93dVfWgqvqVqvrFeRcGAMD8TTWBJMlrk1yZ5JuSpKoeXlWXz7wqAAAWYtoweJfu/s0ktyXJzq4kD595VQAALMS0YfCWqnpgkk6Sqqok95x5VcDy2JUEYFCmWlomyQuS/HKSL6+qpye5NIk9iwEA1tRUYbC7/6iqLk3ypCSPTHImyWXzKAxYsN1u4PZ2cvasNQgBBmLaRaf/RXe/NKNJJFfuHLvrPAoDlmi/XUmEQ4CNNO0w8f2q6nu7+9eSpKouSHJFksfOvDJgsQ7alUTHEGCjTRsGfyjJ26vqpowmkbwqyYtmXhWwPOcKefYxBtg4k25Hd3mS65K8N8lzk/xaks8neVJ33zS/8oCls48xwEabtDP4mowmjDxj5+8Lk7wnyfdV1Q3dfeV8ygNWhhAIsJEmCoPd/Y4k79i9XlXnJbk4o2D4mOxMJgEAYL1Me85gkqS7P5/k+p0/r51pRQAALMy0O5AAALBBhEEAgAETBgEABmzqMFhV1hUEANgQh04gqaorxq8mOZnkp+ZWEQAACzPJbOJPdPczd69U1cvnWA8AAAs0yTDxv9xz/UfmUQgAAIt3aBjs7v+WJFV1/s71j827KGCNbG3dsVUdAGtnmgkkl82tCgAAlmKaHUhqblUA62e3G7i9nZw9O7pu/2KAtTNNZ7DnVgWw3k6cWHYFAByRziBwNLtdwN0Ooa4gwFqaJgy+ZG5VAOtLCARYaxMPE3f3DVV1pqq+NEmq6jlV9YKquvv8ypufqnpSVb2yqt5cVY9fdj0AAMsw7XZ09+7uT1TV1yZ5VpK/nuSVsy/r3Krqsqr6cFXdsOf4pVX1gaq6qapefK7n6O7f6O5nJfkHSf7eHMsFAFhZ0wwTJ8nnquq8JN+f5Ke6+4qqumYOdR3m1Ul+Mcnluweq6q5JXpbkcUluTvKeqnpLkrsm+ck9j39Gd3945/JLdx4HADA404bBn0/yviT3SLLbebvXTCuaQHe/q6ou3HP40Ulu6u4PJUlVvT7JE7v7J5N8x97nqKpK8q+S/GZ3X3fQa1XVs5M8O0ke8IAHzKR+AIBVMdUwcXdfnuTrkzy8uz9dVRclefdcKpve/ZL86dj1m3eOHeT5Sb41yd+tquccdKfufkV3n+ruUxdccMFsKgUAWBHTdgbT3Z8cu3xTkqfPtKKj22/pmwPXRuzun8+o0wnMmuVmANbGtBNIVtnNSb5y7Pr9k9yypFoAANbC1J3BFfaeJA+qqgcm+bMkT03yvcstCQZmvy3qEh1CgBW2lp3BqnpdRucqPriqbq6qH+zuzyd5XpK3J7kxyRXd/QfLrBMGzRZ1AGvh0M5gVT0uyVOSvKy7t6vq2d39ivmXdrDu/p4Djl+d5OoFlwPsskUdwNqZZJj4hzOaJPLSqrpPkpPzLQlYe0IgwNqYZJj4I9398e5+YZLHJ/m6OdcEAMCCTBIG37Z7obtfnLFdPwAAWG+HhsHufnOSVNX5O9d/Yd5FAQCwGNPMJr5sblUAALAU04TB/Xb4AJjM1tYds4wBWBnTLDp94NZuAOe0tWUhaoAVpTMILI6FqAFWzjSdwZfMrQpgs50+rSMIsKIm7gx29w1JUlV3m185AAAs0jSdwVTVLyd5clV9KsktSa5Pcr3lZoBD6QgCrKSpwmCSb0zyZd39uaq6X5JHJvnq2ZcFAMAiTBsG/3OSv57kw939Z0n+LMnVM68KAICFmGY2cZK8IsmZqnphVX1jVZkaCACwxqYNg7+a5IqMOoo/nOQ/VdUHZ14VAAALMe0w8c3d/aPjB6rqi2ZYDzAklpsBWLppO4PbVfWPxg9092dmWA8AAAs0bWfwy5J8a1W9KMl1Sd6XZLu73zDzyoDNZos6gJUwVRjs7qcktw8NPyzJI5J8fRJhEDgaW9QBLNW0i06fSnJjd38qo87gdXOpCth8tqgDWAnTDhNfnuRRu1eq6vwkj+nut860KmAY9guBAiLAQk07geTW7r5190p3fzTJj822JGCwds8jPHPmjlAIwFxNGwY/VFXfvufY3WdVDEAS5xECLNC0w8TPT/KbVfX3M9qa7mFJLDoNzIbzCAEWbqrOYHf/eZKvTfLGJBckuT7J986hLgAAFmDazmC6+7aMwuAbZ18OMHg6ggALNe05gwAAbBBhEFh9W1t3nl289zoARzZVGKyqM1X1pTuXn1NVL6gqs4kBANbUtOcM3ru7P1FVX5vkWUnemuSVSX5g5pUB7Hb/9u5hPH7dOYYAxzJtGPxcVZ2X5PuT/FR3X1FV18yhLoA727v24N7rlqQBOJJpw+DPJ3lfknskefHOsXvNtCKAXbvBbm/QE/wAZmaqMNjdl1fVrye5rbs/XVUXJXn3fEoD2HGu0GfoGOBYjrLO4CfHLt+U5OkzrQjgMAeFPdvYAUxtqjBYVd+c5GlJPp7khox2ILmhuz8zh9oADmfoGOBYpu0M/mqS5+487quTPCmj/YkvmnFdAAAswLRh8KbuftPO5TfMuhiAI9MRBDiSaXcgOVNV/7iqai7VAACwUNN2Bh+W5OFJXlRV1ybZTrLd3bqEAABraNqlZZ6cJFV1z9wRDB+TDR8yrqonJHnCRRc5NRIA2CzV3ZPfueo+Sf5xkvsm+cMkl3f3/5hTbSvn1KlTfc01NlyBtWB2MTBwVXVtd5867H7TnjP4+iR/meSqJH8tye9U1aOPUB8AACtg2nMGv6K7f3rn8lur6t8m+bWMhooBAFgz03YGP1ZVX717pbs/lFGHEGB1bG2Ntqc7c+aO4WIA9jVtZ/CHklxZVf8xye8nuTjJB2deFcAs2J4O4FDThsFPJ3lURjuPXJzkfUleOOuiAI7l9GkTSAAmNG0YfFN3PyrJFbsHquoxSf7zTKsCmDXhEGBfE4XBqnpKRh3BL6mqhyb5r919287Nr8hon2KA1SH0AUxk0s7g7ya5R5JnJvnZJA+uqo8nuSWjoWOA1bU7oeTs2dFlQRHgdhOFwe7+sySXV9UHu/t3k9sXoH5gkvfPsT6A2TGhBOALTHvO4I1V9eO58w4kn5p9WQAztN+EEucQAiSxAwkAwKDZgQQYhvEOoHMIAW5nBxJguJxDCHDsHUgeFjuQAOvGotQAt5uqM9jd789ovcF3ZjSJZDvJ98yhLgAAFmCqzmBVfXOSpyX5eJIbklyfpOdQFwCwrnTe18q0w8S/muS5O4/76oz2KH5YkotmXBfAfPmQAkgyfRi8qbvftHP5DbMuBgCAxZronMGquryqXpDk3VX1f825JgAAFmTSCSSv2bnvlyf5+1X1x1X1lqr68ar67vmVBwDAPE26N/E7krxj93pVnZfk4iSPTPLoGDIGAFhLE4XBqvr+JD+bUXfwqiTP6+7rM5pNDLD+zH4EBmrSYeJ/luRxSR6S5E+S/MTcKgJYRVtbdwRG4GC72z2eOXPHvxn/flbapGHwE9393u7+cHf/04yGhgE2w34fXsDx2O5xbUy6tMxXVNWzk9yY5P1J7ja/kgCW5KAPr92wePbs6LKhZDiY7R7XzqRh8EczWmT6aUkekeReVXV1kvclub67Xzen+gDmb78Pr/0+zHQ6gA006WziV4xfr6r7ZxQOH5Hkf08iDAKbS6cD2GCTziau7r59D+LuvjnJzUmuPug+AGtlPOQZFgYGZNIJJO+squdX1QPGD1bV3avqm6vqNUl+YPblASzRuc4hNNEE2BCTnjN4aZJnJHldVT0wyceT3DOjMPnvkvw/3b09nxIBFmy/YWHdQWBDTXrO4K1JfinJL1XV3ZKcn+TT3f3xeRYHAKy5vaddJP5ztWImHSa+XXd/rrv/XBAEBsuahDA9s/FX1qGdwap6XJKnJHlZd29X1bP3zi4G2DiHdS58sMHBxv/9mI2/8iYZJv7hJE9P8tKquk+Sk/MtCWDFndx5G/TBBmyASYaJP9LdH+/uFyZ5fJKvm3NNAAAsyCRh8G27F7r7xUkun185AAAs0qHDxN395qp6SHe/f+f6L8y/LIAVdtThYedNASto0tnEV1fVZXsXnQYAYL1NGgYfkuS9Sc5U1c9V1QVzrAkAGCo7/CzcRGGwuz+7Mzz80Iz2JP69qvqxqvqSuVYHsEmsT8hQnT7t9IgVNtWi0919a3f/TJJHJLk1yXVV9cK5VAawiaxPCKyYSfcmTpJU1YUZDRk/OMkDkvxlkp9I8jOzLgxgbUw6McT6hMAKmigMVtX1Se6f5E+SvD/JjUl+O8nLknxgbtUBrJu9wdAMYmDFTdoZ/DtJPtTdPc9iAABYrInCYHd/MEmq6kFJXpzk0939vHkWBrAWtrZGE0POnv3Cy3u7gXuvj3cNdRCBJZlqAkmS1ya5Msk3JUlVPbyq7EgCMD4xZPfybjg87gxiS20wBLu/5+P/ZvzeL8S0YfAu3f2bSW5Lku6+IcnDZ14VwLo4fXo0MeTkyS+8vMsMYpiOfzMLNdVs4iS3VNUDk3SSVFUluefMqwLYBJMO/x42vAxDYNLV0kwbBl+Q5JeTfHlVPT3JpUlumHlVAOvkKB9a4wHwkktGx/brhgiKwJxNGwafmVEAfFKSRyY5k+Q1sy5q1VTVE5I84aKLLlp2KcCqmzasnTixfwdx77lShs2AOZk2DN4vyXd39+uSXFlV903y1iSPnXllK6S7r0py1alTp5617FqANXNQOBwPgIc93rAZMEfThsEfSvL2qvpgkr9K8uokL5p1UQAALMakO5BcnuS6JO9N8twkv5bk80me1N03za88gIEY7/qNnycIm0zXeyVM2hl8TUbnCD5j5+8Lk7wnyfdV1Q3dfeV8ygPYYId9AE5ynuBxPkx9EAOZfAeSdyR5x+71qjovycUZBcPHZLQQNQCzsPd8QmENmKNpzxlMknT355Ncv/PntTOtCAAYLv/5WbhpdyABAGCDHKkzCMAKONdEk4POB3SeILCHMAiwig4Ka/uFuaMuSG1nEyDCIMD6mnTh6nOxswnLsPt76z8kK0EYBBiqkydHf/sQZpn8h2TphEGATaTjwirbuwe339GlEgYB1sm5Qt7eD9b9Oi7jjweIMAiwfsZD3kEdlXMNARuWA8YIgwDr5Ljn+Y0/XncQiDAIMFzO0wIiDAJshvEFqHX8gCkIgwDr5LBu3u75gLp+wISEQYBNMH4O4LmCoJDIKvH7uBLusuwCAABYHmEQADi6rS3nqa45YRAAYMCcMwiwKeZx/pXtwmDj6QwCbLp5DOMZGoSNIQwCMHvCIqwNYRAAYMCEQQB08mDATCAB2HTj29RNMhFkmkkjJpgM1+7Pfu82iH4X1o4wCDAEu9vUTWv3gz7xYc/Bjvr7xUoQBgE23cmTo7+PGuL2ftBvbd0REi+55DiVsc52f5/8J2HtCYMAHGw3SO7noG7QeFgEVp4wCMB0Tp+ebLKJoUNYC8IgAHeeBLB3QsC0Jg2LrBfDwRvL0jIAjIx38ibt6m1vJ2fOCH+wxnQGATbdJJ2c8Ukm4x0gw8Gw8YRBAA530KSQ485UZv352a89w8QAHOz06Tt/2OsCwsbRGQTgcCaFwMbSGQRgens7hqw3e1MPms4gALMPdoIirA1hEIA7O2qQsw7d5tpv/clN+Dlv0tdyDMIgAEc38A/RwVmHCUQC3tSEQQCWb78P8Hl9qAsL09u7/iQbRRgEYDKzCgF7JyoMPVzsDVlCFwtmNjEA87Mps1Q35euAfegMAsBQ7Qbc8ckhOpKDozMIwGrTlVuMdZgcsoo24PdTZxCA2dhv72JW25DOU9wvsO1dLmevTf5+jNEZBGB25t1dOqwLM8nt29vJmTNr382Zuw3oeE1ER1RnEIAZOXny6I/dr6u4e+ySS45X134EgOmtapdsb2A9qNu3X/1D6IhOQBgEYLEOCnn7BbRZh7ajrJe3d/eNWQeHveHFZI7jE/anIgwCMD8HnUd44sSdQ89uV3E8qB2n03hc+wXGeQaMvc8tzExub2DW7ZuaMAjAfO0Gm90P6d2Qt729HqFn3rtv7H4/hjSZg5UiDAIwP+fq7o2HoHNNVDhsxueuvcO5u8/N5lj0z3WSCTST/n4mK/v7KAwCsPqm6SCuQ7fxOOYRiFY0pKyNNf+dEwYBWG3jHcR73/vgtQyPM5yrk8h+Jvl92IDfHWEQgNlY1Ifhcbswa/yhvRBDmtW8qV/XlITBCVTVE5I84aKLLlp2KQDDtswZxutmFh2rNR/+ZDLC4AS6+6okV506depZy64FYKPt14067i4Ys1hXcJM6SId9P9Z9VvO61bsChEEA5m88YHXfcWzX7gf4ve+9Ot2oc9Uxy6A0i+dY1+A2NCv6cxIGAZifvV2m8YC13wfi3jX3zmVeH6yzWFfwsMef6/Zzvablc5gDYRCA+TtsLcFZvs401nWyxKp0Txdh/PdmvzX91uHnteKEQQBW27rMUl6Eee+Gsg7W4ee0ZoRBADbXfpNBxq3DZIlVrm0Rxr/uVfxerFItRyQMArD5htZN2u/cwl3rNiTO3AmDAKyPacPL0IdVDwrBk4TjIX6/BkoYBGD9zHtW7TQzevfe97DbZ+2g78VB35chh2P2dZdlFwAAR7bM4d/DXvug23fD25kzs51hPbShcGZGZxCA1TFpt2qZw7+HvfZ+t+93/6OGt/0mwej2cQzCIACLscpB5bi1TRPGhDdWjDAIAKwH4XkuhEEANts6B4hFT0ZhkEwgAYBVZ3IIc6QzCMB6GkKX7DjnFx50/yF835iKMAgA607A4xiEQQCY1t79cvee13fY7bBChEEAmIVpFqEWCFkhwiAAHMdRFqGGFWI2MQDAgAmDAAADJgwCAAyYcwYBYJU5z5A50xkEABgwYRAAYMAMEwPAcR02lGuolxWmMwgAMGDCIADAgAmDAAADJgwCAAyYMAgAMGBmEwMALML2dnL27LKr+AI6gwAAi3LixLIr+AI6gwAAi3Dy5LIr2JfOIADAgAmDAAADJgwCAAyYMAgAMGDCIADAgAmDAAADJgwCAAyYMAgAMGDCIADAgAmDAAADJgwCAAyYMAgAMGDCIADAgAmDAAADJgwCAAyYMAgAMGDCIADAgAmDAAADJgwCAAyYMAgAMGDCIADAgAmDAAADJgwCAAyYMAgAMGDCIADAgAmDAAADJgwCAAzY4MNgVT20qv51VV1ZVf9w2fUAACzSwsNgVd17J3i9v6purKq/fcTnuayqPlxVN+xz26VV9YGquqmqXnyu5+nuG7v7OUmekuTUUWoBAFhXy+gM/r9Jfqu7H5LkkUluHL+xqu5bVV+y59hF+zzPq5NcuvdgVd01ycuSfHuSi5N8T1VdXFWPqKq37vlz353HfGeS30nyjuN/eQAA62OhYbCqvjTJNyX5lSTp7s9298f33O2SJG+uqnvsPOZZSX5+73N197uSfGyfl3l0kpu6+0Pd/dkkr0/yxO7+/e7+jj1/PrzzXG/p7v81ydNm9KUCAKyF8xb8el+V5CNJXlVVj0xybZJ/1N2f2r1Dd7+hqh6Y5PVV9YYkz0jyuCle435J/nTs+s1Jvv6gO1fVVpInJ/miJFcfcJ8nJHnCRRft16AEAJjA6dPLrmBfix4mPi/Jo5K8vLu/JsmnknzBOX3d/dNJbk3y8iTf2d2fnOI1ap9jfdCdu/t0d/+f3f1D3f2yA+5zVXc/+8SJE1OUAQCw+hYdBm9OcnN3/97O9SszCod3UlXfmOThSd6U5EeP8BpfOXb9/klumb5UAIDNt9Aw2N1/keRPq+rBO4e+Jckfjt+nqr4mySuTPDHJ05Pcp6r+xRQv854kD6qqB1bV3ZM8Nclbjl08AMAGWsZs4ucn+TdVdX2Sk0l+Ys/tfy3Jd3f3B7v7r5L8QJI/3vskVfW6JO9O8uCqurmqfjBJuvvzSZ6X5O0ZzVS+orv/YG5fDQDAGqvuA0+nY49Tp071Nddcs+wyAAAOVVXXdvehaygPfgcSAIAhEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAasunvZNayNqvpIkj+e8O4nkpydYzmztAq1LrKGeb7WrJ97Fs933Oc4P8lHj1kD87EK/3YXaZ2+3lWoddE1zOv1VvF9dRbPs4j31r/Z3RccdidhcE6q6hXd/exl1zGJVah1kTXM87Vm/dyzeL7jPkdVXdPdp45TA/OxCv92F2mdvt5VqHXRNczr9VbxfXUWz7NK762GiefnqmUXMIVVqHWRNczztWb93LN4vlX4+TIfQ/vZrtPXuwq1LrqGeb3eKr6vzvJ5lk5nEDinVfrfK8CmWKX3Vp1B4DCvWHYBABtoZd5bdQYBAAZMZxAAYMCEQQCAARMGAQAGTBgEJlZVX1xVr6mqV1bV05ZdD8AmqKqvqqpfqaorl/H6wiAMXFVdVlUfrqob9hy/tKo65dhoAAAE70lEQVQ+UFU3VdWLdw4/OcmV3f2sJN+58GIB1sQ0763d/aHu/sHlVCoMAsmrk1w6fqCq7prkZUm+PcnFSb6nqi5Ocv8kf7pzt9sWWCPAunl1Jn9vXSphEAauu9+V5GN7Dj86yU07/1v9bJLXJ3likpszCoSJ9w+AA0353rpU3syB/dwvd3QAk1EIvF+SX0/yXVX18mzQVkwAC7Lve2tV/Y2q+tdJvqaqXrLoos5b9AsCa6H2Odbd/akkT190MQAb4qD31v+e5DmLLmaXziCwn5uTfOXY9fsnuWVJtQBsipV8bxUGgf28J8mDquqBVXX3JE9N8pYl1wSw7lbyvVUYhIGrqtcleXeSB1fVzVX1g939+STPS/L2JDcmuaK7/2CZdQKsk3V6b63uXnYNAAAsic4gAMCACYMAAAMmDAIADJgwCAAwYMIgAMCACYMAAAMmDAKDV1W3VdX22J8XL7umXVV1ZVV91c7lT07xuNNV9W17jr2gqn6pqi6oqt+ada3AerI3MUDy6e4+OcsnrKrzdhaYPc5zPCzJXbv7Q0d4+Osy2t3g7WPHnprkn3T3R6rqz6vqG7r7d49TI7D+dAYBDlBVf1RV/7yqrquq36+qh+wc/+Kquqyq3lNV762qJ+4c/wdV9YaquirJv6uqu+x04v6gqt5aVVdX1d+tqm+pqjeNvc7jqurX9ynhaUnevE9d51fVu6vq/9i5/k92arm+qv75zt2uTPIdVfVFO/e5MMn/kuR3dm7/jZ3nBwZOGARI7rlnmPjvjd320e5+VJKXJ3nhzrEfSfLb3f11SR6b5P+uqi/eue1vJ/mB7v7mJE9OcmGSRyR55s5tSfLbSR5aVRfsXH96klftU9c3JLl2/EBVfVmStyX5Z939tqp6fJIHJXl0kpNJvraqvqm7/3uS/5Lk0p2HPjXJv+07tp26Jsk3Tvj9ATaYYWKAcw8T73bsrs0o3CXJ45N8Z1XthsN7JHnAzuV/390f27n8vyV5Q3f/VZK/qKp3Jkl3d1W9Nsn3VdWrMgqJ37/Pa39Fko+MXb9bknckeW53nxmr5fFJ3rtz/V4ZhcN35Y6h4jfv/P2Msef6cEadQmDghEGAc/vMzt+35Y73zEryXd39gfE7VtXXJ/nU+KFzPO+rklyV5NaMAuN+5xd+OqOguevzGYXSb0uyGwYryU929/+3z+N/I8nPVtWjktyzu68bu+0eO88PDJxhYoDpvT3J86uqkqSqvuaA+/1Oku/aOXfwy5Js7d7Q3bckuSXJS5O8+oDH35jkorHrnVF37yFjM57fnuQZVXWvnVruV1X33XmNTyY5neSyjLqE4/5WkhsO+0KBzaczCLBzzuDY9d/q7nMtL/PjSX4uyfU7gfCPknzHPvd7Y5JvySh0/dckv5fk7Njt/ybJBd39hwe8ztsyCpD/YfdAd99WVU9NclVVfaK7f6mqHprk3TvZ9JNJvi+jYeBkFAJ/PaNh4nGP3Xl+YODqjnOJAZi1qrpXd3+yqv5GRhM6vqG7/2Lntl9M8t7u/pUDHnvPJO/cecxtM67rXUme2N3/Y5bPC6wfYRBgjqrqdJJ7J7l7kp/u7lfvHL82o/MLH9fdnznH478tyY3d/SczrOmCjALmb8zqOYH1JQwCAAyYCSQAAAMmDAIADJgwCAAwYMIgAMCACYMAAAMmDAIADNj/BPgNX5a+SacjAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 720x576 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "data06 = pd.read_csv(\"energy_2006.csv\")  \n",
    "x06 = data06['x']\n",
    "y06 = data06['y']\n",
    "xerr06 = data06['xerr']\n",
    "yerr06 = data06['yerr']\n",
    "\n",
    "\n",
    "plt.figure(figsize=(10,8))\n",
    "plt.loglog()\n",
    "plt.errorbar(x06, y06, yerr=yerr06, xerr=xerr06, fmt=',r')\n",
    "plt.xlabel(\"Energy (keV)\")\n",
    "plt.ylabel(r\"$keV^2 (Photons{\\ }cm^{-2}{\\ } s^{-1}{\\ } keV^{-1})$\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoMAAAHnCAYAAADD41dfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3X2UrWdZH+DfTQKCYo/SBLEBmtRQISQQMCKtxR6pYLR8iUpBEBsERMFVrFjC0kpRq6KttSpig8QQrKQhgCQQTTVyQrGIJOEQg4GukArG4AqIHD4WIIS7f+w9ZDKdObP3nP0573WttRfzvrPn3feQdfb89v28z/NUdwcAgGG607ILAABgeYRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwI5fdgHr5IQTTuiTTz552WUAAOzqmmuu+Uh3n7jb84TBKZx88sm5+uqrl10GAMCuquoDkzzPMDEAwIAJgwAAAyYMAgAMmDAIADBgwiAAwIAJgwAAAyYMAgAMmDAIADBgwiAAwIAJgwAAAyYMAgAMmDAIADBgwiAAwIAJgwAAAyYMAgAMmDAIADBgwiAAwIAJgwAAAyYMAgDMy8GDo8cKEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwI5fdgEAAPvKwYO3f334cHLkyB3PHTq04IKOTmcQAGCeDhxYdgVHpTMIADBLmzt/Gx3BFesGbqYzCAAwYMIgAMCACYMAAAMmDAIADJgwCAAwYMIgAMCACYMAAAMmDAIADJgwCAAwYMIgAMCACYMAAAMmDAIADJgwCAAwYMIgAMCACYMAAAMmDAIADNjgw2BVPaCqfqOqLqmqH1x2PQAAi7TQMFhVd62qP62qd1fVe6rqJcdwrfOr6taqun6b751dVe+rqhur6tyjXae7b+ju5yR5UpKz9loPAMA6WnRn8LNJHtndD05yZpKzq+rhm59QVfesqi/fcu7Uba51QZKzt56squOSvCzJtyU5LclTquq0qjqjqt605XHP8c88Lsnbklx57L8iAMD6WGgY7JFPjg/vPH70lqf98yRvrKq7JklVPSvJr2xzrbcm+eg2L/OwJDd2903d/XdJLkry+O7+s+5+zJbHreNrXdrd/zTJU7eru6oeW1XnHTlyZPpfGgBghS38nsGqOq6qDie5NckfdPc7Nn+/u1+b5PeTXFRVT03yjIyGcCd1UpK/3HR88/jcTvUcrKpfqar/luTy7Z7T3Zd197MPHDgwRRkAAKvv+EW/YHffluTMqvqKJG+oqtO7+/otz/mFqrooycuTfM2mbuIkaruXPUo9h5IcmuL6AAD7xtJmE3f3xzIKYdvd9/eIJKcneUOSF0956ZuT3GfT8b2T3LK3KgEA9rdFzyY+cdwRTFXdLcm3JHnvluc8JMkrkjw+yTlJ7lFVPzPFy7wzyf2q6pSqukuSJye5dBb1AwDsN4vuDH51krdU1XUZhbY/6O43bXnOlyb57u5+f3d/Icn3JfnA1gtV1WuSvD3J11bVzVX1/UnS3Z9P8rwkVyS5IcnF3f2euf1GAABrrLp3vJ2OLc4666y++uqrl10GAMCuquqa7t51DeXB70ACAByDgwdHD9aWMAgAMGDCIADAgAmDAAADJgwCAAzYwncgAQBW3DQTQg4fTo4cmfxnDh3aQ0HMk84gAHBsDhxYdgUcA51BAOCOpunebXQEdfzWls4gAMCACYMAAAMmDAIADJh7BgGAvXOv4NrTGQQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAVt/Bg6MHMycMAgAMmDAIADBgwiAAwIAJgwAAAyYMAgAM2PHLLgAAGJC9zgg+fDg5cmRvP3/o0N5ecyB0BgGA9XDgwLIr2Jd0BgGAxdlrl26jI6jLN3M6gwAAAyYMAgAMmDAIADBgwiAAwIAJgwAAAyYMAgAMmKVlAIDVZ0mZudEZBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYsD2Fwar6sqo6btbFAACwWBOFwaq6U1V9T1W9uapuTfLeJB+qqvdU1S9W1f3mWyYAsCcHD44esINJO4NvSfI1SV6U5F7dfZ/uvmeSRyT5kyQ/X1VPm1ONAACLN5AgffyEz/uW7v7c1pPd/dEkr0vyuqq680wrAwBg7ibqDG4OglX1wt2eAwDAeti1M1hVF28+THJmkpfOrSIA4OimGbo8fDg5cmS6nzl0aMqCWGeTDBN/vLufuXFQVS+fYz0AwKwdOLDsClhhk4TB/7jl+MfnUQgAMKFpOncbHUHdPnaw6z2D3f1/k6SqThgff3TeRQEAsBjTLDp9/tyqAABgKaYJgzW3KgAAWIppwmDPrQoAAJZi0kWnE51BAFg/Jo6wi2k6gy+aWxUAACzFRGGwqu7f3dfPuxgAABZr0mHiy6vqUJL/0N0fnGM9AACzNc3uK5vtZfeWZO2G5icdJr5/kncluaqqfrmqTpxjTQAAq2EAu7dM1Bns7r9L8qtV9Yokz0vyjqr67SS/2N2fmGeBAADHZK+duoHs3jLNBJJ092e6+z8lOSPJZ5JcW1UvmEtlADAkBw/ufTgTjsFUYbCqTq6qs5M8M8l9k3wiyc/OozAAAOZvomHiqrouyb2TfDDJe5PckOSPkrwsyfvmVh0AAHM16Wzi70hyU3fbhQQAYB+ZaJi4u9/f3V1V96uqV1bVr827MAAA5m+a7eiS5NVJXpLkpUlSVacn+Xfd/fRZFwYA+8Kkk0KmWdNun89uZbGmmkCS5E7d/XtJbkuS8a4kp8+8KgAYogGsacfqmbYzeEtVnZKkk6SqKsndZl4VAOwXk3bxBrKmHatn2jD4/CS/meReVXVOkrOT2LMYAGBNTRUGu/svxusMPiHJg5NcleT8eRQGAMD8TRUGq+pnuvsnklwyfqSqjptHYQAAzN+0E0hOqqrv2TioqhOT/OFsSwIAYFGmvWfwB5JcUVU3ZjSJ5LeSvHDmVQHA0Jg4wpJMuh3dhUmuTfKuJM9N8jtJPp/kCd194/zKAwBYkoEE9EmHiV81fu4zMgqCJyf52yRPq6rvmk9pAADM20Sdwe6+MsmVG8dVdXyS0zKaUfzwjCeTAACwXqa9ZzBJ0t2fT3Ld+PHqmVYEAKvOAtHsI9POJgYAYB8RBlkvBw9Ovuk7ALCrwYfBqnpAVf1GVV1SVT+47HoAABZp6nsGq+qF3f3SvbxYVd0nyYVJ7pXkC0nO6+7/usdrnZ/kMUlu7e7Tt3zv7CT/NclxSX6zu39+p+t09w1JnlNVd0ryir3Uwhxt7QIePpwcObJ9d9C9OwAwtV3DYFVdvPkwyZlJ9hQGM1qb8Ee7+9qq+vIk11TVH3T3n296vXsm+XR3f2LTuVO3Wc/wgiS/llG43FzvcUleluRRSW5O8s6qujSjYPhzW67xjO6+taoel+Tc8fVYdQcOLLsCYD+b5FaUo30w3coHVVbcJJ3Bj3f3MzcOqurle32x7v5Qkg+Nv/5EVd2Q5KQkf77paf88yQ9W1bd392eq6llJviPJt2+51lur6uRtXuZhSW7s7pvG9V6U5PHd/XMZdRK3q+vSJJdW1ZszWkeRVbH1TdQMPmBV+GDKPjFJGPyPW45/fBYvPA5yD0nyjs3nu/u1VXVKkouq6rUZLXT9qCkufVKSv9x0fHOSbzhKHQeTPDHJlyS5fIfnPDbJY0899dQpygBgLU3yYdMHU/aRXcNgd//fJKmqE7r7I9390WN90aq6e5LXJXl+d398m9f8hXFH7+VJvqa7PznN5bc51zs9ubsPJTl0tAt292VJLjvrrLOeNUUdzIM3XgCYqWlmE58/ixesqjtnFAT/e3e/fofnPCLJ6UnekOTFU77EzUnus+n43klu2UOpAAD73jRhcLuO21SqqpK8MskN3f1LOzznIRnN6n18knOS3KOqfmaKl3lnkvtV1SlVdZckT05y6bFVDgCwP00TBnccap3CNyb53iSPrKrD48e3b3nOlyb57u5+f3d/Icn3JfnA1gtV1WuSvD3J11bVzVX1/ckXt8p7XpIrktyQ5OLufs8MageAkUOH3LbCvjHNOoPH3Bns7rftdp3u/uMtx5/LNuv/dfdTjnKNy7PDZBAAAG43TWfwRXOrAgCApZg4DHb39VV1VVX9vSSpqudU1fPH9+UBALCGpt2b+Cu6++NV9XVJnpXkK2MLNwCAtTXt3sSfq6rjkzw9yUu7++KqunoOdQEAsADThsFfSfLuJHfNaC/fJLn7TCsCAGBhpgqD3X1hVb0+yW3d/emqOjWj5V0AAFhD03YGs3lruO6+MaOFoQEAWEPTTiABAGAfEQYBAAZMGGQ1HDw4egAAC7VrGKyqR1XVK6rqzPHxs+dfFgAAizDJBJIfymiSyE9U1T2SnDnfktj3tusAHj6cHDny/3/PRvAAMFeTDBN/uLs/1t0vSPLoJF8/55oYqgMHZnctw84AMJFJOoNv3viiu8+tqh+eYz0MwXbdvo3gphMIzIP3GNjRrmGwu9+YJFV1Qnd/pLt/df5lwR5s7gQadgaAiUyz6PT5SR43r0IYuHmEtFkOOwPAPjVNGKy5VQGzsDlQGhICgIlMs85gz60KAACWQmeQ/UlHEIZjkpUDdrqPeDvePxiYaTqDL5pbFQAwb+4jhm1N3Bns7uuTpKru3N2fm19JADCFSTp57iOGHU0zTJyq+s0kT6yqTyW5Jcl1Sa6z3AwAwHqaKgwmeUSSr+ruz1XVSUkenORBsy8LAIBFmDYM/kmSr0xya3f/VZK/SnL5zKsCAGAhpg2D5yW5qqpemeQdGQ0RH5l9WQAwQ+4VhB1NM5s4SX47ycUZhcgfSvK/q+r9M68KAJLRxI9JloMB9mzazuDN3f3izSeq6ktmWA8AAAs0bWfwcFX9m80nuvuzM6wHAIAFmrYz+FVJvqWqXpjk2iTvTnK4u18788oAAJi7qcJgdz8p+eLQ8AOTnJHkG5IIg9yRBV6BSex2P+Ak28h5n4FjMu2i02cluaG7P5VRZ/DauVQFABtsIwdzNe0w8YVJHrpxUFUnJHl4d79pplWxXrb7xH60T/M+xQMbdns/MMoAczftBJLPdPdnNg66+yNJfmq2JbFv7IdP85a1AGCfm7YzeFNVfVt3/96mc3eZZUGsoe0+sfs0DwBrYdow+MNJfq+qvjejrekemMSi0wAAa2ra2cQfqqqvS/KEJGcmuS7Jv51HYay5de4Ibh4W3unex3X+/QBgk2k7g+nu25K8bvyA/W8/3PsI68oHL5i7qcMg7Hub//i49xGAfW7a2cSwv5gtDMDATRUGq+qqqvp746+fU1XPryqziQEA1tS0ncGv6O6PjyeRPCvJVyZ5xezLghVx6JAhYgD2tWnvGfxcVR2f5OlJXtrdF1fV1XOoC+Zj65DwdrOFhT8ABmTaMPgrSd6d5K5Jzh2fu/tMK4JFM1sYgAGbdp3BC6vq9Ulu6+5PV9WpSd4+n9JgDrZ2/WYxW9iMYwDW2F7WGfzkpq9vTHLOTCsCAGBhpgqDVfXIJE9N8rEk12e0A8n13f3ZOdQG87fXbt7RdinRIQRgjUzbGfztJM8d/9yDMtqW7oFJTp1xXbBe3HcIwJqaNgze2N1vGH/92lkXA2vDLiUA7BPTrjN4VVX9SFXVXKphNdmlAwD2rWk7gw9McnqSF1bVNUkOJznc3bqEAABraNqlZZ6YJFV1t9weDB8eQ8b7w07dv+0WZt5s6MOjQ//9AVhr084mvkeSH0lyzyR/nuTC7r5gDnWxakyQAIB9adph4ouS/GGSdyQ5I8nbquqc7v7TmVfG4u3U4TJBAgD2rWnD4Fd39y+Mv35TVf2PJL+T0VAxAABrZtrZxB+tqgdtHHT3TUm+dLYlAQCwKNN2Bn8gySVV9b+S/FmS05K8f+ZVsVoMDwPAvjVtZ/DTSR6a5C0ZTSJ5d5KnzLooAAAWY9rO4Bu6+6FJLt44UVUPT/InM60KAICFmKgzWFVPqqqfT/LlVfWAqjpu07fPm09pAADM26SdwT9Octckz0zyS0m+tqo+luSWjIaOAQBYQxOFwe7+qyQXVtX7u/uPky8uQH1KkvfOsT4AAOZo2nsGb6iqn84ddyD51OzLgn3Cgt0ArLhpZxNflOQTSS7LaH3Bt1XVw2ZeFQAAC2EHEgCAAZs2DH60qh7U3dclox1IqsoOJLDZxtBwkhw+nBw5csdzhowBWCHHugPJA2MHEji6AweWXQEsnvtlYW1MFQa7+71V9dAkT8hoK7rDSX50HoXB2tr8x88fRABW3FRhsKoemeSpST6W5Pok1yXpOdQFAMACTDtM/NtJnjv+uQdl1CF8YJJTZ1wXAAALMG0YvLG73zD++rWzLgb2HcPDAKy4icJgVV2Y5Nokb6+qH+3u/zzfsgBYaZtnyG9nu5n02/GBCZZu0kWnXzV+7r2SfG9VfaCqLq2qn66q755feQCsLTPpYS1MujfxlUmu3DiuquMzmk384CQPiyFjgGHZraNnJj2sjUmHiZ+e5Jcy6g5eluR544Wnr5tjbQAAzNmkw8Q/meRRSe6f5INJfnZuFQEAsDCTzib+eHe/a/z1v6+qd8yrIAAAFmfSMPjVVfXsJDckeW+SO8+vJAAAFmXSMPjijBaZfmqSM5LcvaouT/LuJNd192vmVB8Ay3CsE0BMHIG1Mels4vM2H1fVvTMKh2ck+fYkwiAAwBqadDZxdfcX9yDu7puT3Jzk8p2eAwDA6pt0NvFbquqHq+q+m09W1V2q6pFV9aok3zf78gAAmKdJ7xk8O8kzkrymqk5J8rEkd8soTP7PJP+luw/Pp0QAAOZl0nsGP5Pk15P8elXdOckJST7d3R+bZ3EAzNHR9g3ebW9hE0Rg35i0M/hF3f25JB+aQy0ArBJ7C8Mg7BoGq+pRSZ6U5GXdfbiqnr11djEAa+ho3T17C8NgTNIZ/KEk5yT5iaq6R5Iz51sSAACLMsls4g9398e6+wVJHp3k6+dcEwAACzJJGHzzxhfdfW6SC+dXDgAAi7TrMHF3v7Gq7t/d7x0f/+r8ywJgqdwrCIMx6aLTl1fV+VsXnQYAYL1NGgbvn+RdSa6qql+uqhPnWBMAAAsyURjs7r8bDw8/IKM9id9RVT9VVV8+1+oAAJirSTuDSUY7kXT3f0pyRpLPJLm2ql4wl8oAAJi7qcJgVZ1cVWcneWaS+yb5RJKfnUdhzMHBg0fffgoAGJyJtqOrquuS3DvJB5O8N8kNSf4oycuSvG9u1QEAMFeT7k38HUlu6u6eZzEAACzWRGGwu9+fJFV1vyTnJvl0dz9vnoVxDHYaCj58ODlyZOfvW1cMAAZnqnsGk7w6ySVJvilJqur0qrIjyTo5cGDZFQAAK2TSYeINd+ru36uqn02S7r6+qk6fQ10ci506fBsdQR1AAGBs2s7gLVV1SpJOkqqqJHebeVUAACzEtJ3B5yf5zST3qqpzkpyd5PqZV8V86AjC8BgRAHYxbRh8ZkYB8AlJHpzkqiSvmnVRAAAsxrRh8KQk393dr0lySVXdM8mbknzzzCsDAGDupg2DP5Dkiqp6f5IvJLkgyQtnXRQAAIsx6Q4kFya5Nsm7kjw3ye8k+XySJ3T3jfMrD4CJ7GV9UfcRApl8NvGrxs99RkZB8OQkf5vkaVX1XfMpDYCZsL4ocBST7kByZZIrN46r6vgkp2U0ieThGS1EDcCyWF8U2KNp7xlMknT355NcN368eqYVAQCwMNMuOg0AwD4iDAKsk4MHd54sArAHexomBmBNuFcQ2IXOIADAgAmDsCyG+wBYAcIgAMCAuWcQYJVt7R7vtKOIewOBPRIGYdE2/ohv90fdH3QmYUcRYIaEQVgmf9TZzdYPCHYUAWZMGIRF2/gj7o86ACvABBIAgAETBgEABswwMcA6cVsBMGPCICyLP+oArADDxACrwI40wJIIgwAAAyYMAgAMmDAIADBgJpAALMvmewTtOQwsic4gwKqwPSGwBDqDAMuyuetne0JgSXQGAQAGTBgEABgwYRAAYMDcMwiwCtwrCCyJziAAwIAJgwAAAyYMAgAMmDAIADBgwiAAwIAJgwAAAyYMAgAMmDAIADBgwiDAPB08OHoArChhEABgwIRBAIABEwYB9sLwL7BPHL/sAgD2la0B8fDh5MiR7YPjoUMLKAjg6HQGAebtwIFlVwCwI51BgFna2u3b6AjqAgIrShgE2M12Q7w7Df8KfcCaMUwMsFeGf4F9QGcQYDfbdfsM/wL7hDAIME/CIrDiDBMDAAyYMAgAMGCGiQH2wvAvsE/oDAIADJgwCAAwYMIgAMCACYMAAAMmDAIADJgwCLDZwYPb70UMsE8JgwAAAyYMAgAMmDAIADBgdiAB2HyP4OHDyZEjdzxntxFgH9MZBNjqwIFlVwCwMDqDAJs7fxsdQd1AYCB0BgEABkwYhFVkrTsAFkQYBIZF0Aa4A2EQAGDATCAB2MzEEWBghEFYJRvDl9a6A2BBhEFYVda6mx2LSgPsSBiEVbIRSqx1N1+CNsAXCYPA/mdRaYAdmU0M7B+WjQGYmjAIADBgholhFRnCBGBBdAYBAAZMZxAYFl1XgDsQBteNmZBwR7utIZj49wJwFIaJgf3FGoIAU9EZXDW7LYuxU+djM10QhsQaggDHRGdwHel8AAAzojO4anbraOh8AAAzpDMIADBgOoPrRkcQdubfB8DUdAYBAAZssGGwqh5QVb9RVZdU1Q8uux4AgGVYyzBYVedX1a1Vdf2W82dX1fuq6saqOvdo1+juG7r7OUmelOSsedYLALCq1jIMJrkgydmbT1TVcUleluTbkpyW5ClVdVpVnVFVb9ryuOf4Zx6X5G1Jrlxs+QAAq2EtJ5B091ur6uQtpx+W5MbuvilJquqiJI/v7p9L8pgdrnNpkkur6s1Jfmd+FcMMWV4IgBlayzC4g5OS/OWm45uTfMNOT66qg0memORLklx+lOc9O8mzk+S+973vLOqE+duPgXE//k4AK2A/hcHa5lzv9OTuPpTk0G4X7e7zkpyXJGedddaO1wOOgaAHsDTres/gdm5Ocp9Nx/dOcsuSagEAWAv7qTP4ziT3q6pTkvxVkicn+Z7llgQztNE9O3w4OXLk9uPtHO05um8AbLKWYbCqXpPkYJITqurmJC/u7ldW1fOSXJHkuCTnd/d7llgmzM+BA7N5zqrbHGZ3CrjCLcAxWcsw2N1P2eH85TnKZBBYaxuhZ5L761b9HrytgW6SbmeyPwIuwIpZyzAI7EPbBb3NYXbVAy7AmhIGgcXbGugEPYClEQZhP5o2VAljAIO1n5aWAQBgSjqDsG6G2r0b6u8NMGfCIAzNdjN2l71si6AHsDSGiYERy7YADJLOIAzNdl04E0gABktnEABgwIRBAIABEwYBAAbMPYOAewUBBkxnEJiPgwe3X8YGgJUiDAIADJgwCAAwYMIgAMCAmUACzMbW+wN32uIuMWEFYIXoDALzY4s7gJWnMwhMZ6et67Ye2+IOYC3oDAIADJgwCAAwYIaJgfkwPAywFoRBYGfbzQTeaZaw8AewlgwTA9MzSxhg39AZnEBVPTbJY0899dRllwKLtV23zyxhgH1FZ3AC3X1Zdz/7gG4IALDPCIMAAANmmBiYjuFhgH1FZxAAYMCEQQCAARMGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGrLp72TWsjar6cJIPzOnyB5IcmdO1Z2nZdS7y9ef1WrO+7iyudyzXOCHJR47x9ZmtZf87XaR1+V2XXaf3zvlcb9XfO/9hd5+425OEwRVRVed197OXXcdull3nIl9/Xq816+vO4nrHco2qurq7zzqW12e2lv3vdJHW5Xdddp3eO+dzvf3y3mmYeHVctuwCJrTsOhf5+vN6rVlfdxbXW/Z/V2ZrSP891+V3XXad3jvnc71l/3edCZ1B4Jis0qdbgHWxSu+dOoPAsTpv2QUArKGVee/UGQQAGDCdQQCAARMGAQAGTBgEABgwYRCYqar6sqp6VVW9oqqeuux6ANZBVf2jqnplVV2y6NcWBoFdVdX5VXVrVV2/5fzZVfW+qrqxqs4dn35ikku6+1lJHrfwYgFWxDTvnd19U3d//zLqFAaBSVyQ5OzNJ6rquCQvS/JtSU5L8pSqOi3JvZP85fhpty2wRoBVc0Emf+9cGmEQ2FV3vzXJR7ecfliSG8efZv8uyUVJHp/k5owCYeI9BhiwKd87l8YbNbBXJ+X2DmAyCoEnJXl9ku+sqpdnn2zVBDBD2753VtXfr6rfSPKQqnrRIgs6fpEvBuwrtc257u5PJTln0cUArImd3jv/JslzFl1MojMI7N3NSe6z6fjeSW5ZUi0A62Ll3juFQWCv3pnkflV1SlXdJcmTk1y65JoAVt3KvXcKg8Cuquo1Sd6e5Gur6uaq+v7u/nyS5yW5IskNSS7u7vcss06AVbIu753V3ct8fQAAlkhnEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwaBwauq26rq8KbHucuuaUNVXVJV/2j89Sen+LlDVfWtW849v6p+vapOrKrfn3WtwHqyNzFA8unuPnOWF6yq48eLyx7LNR6Y5LjuvmkPP/6ajHY2uGLTuScn+bHu/nBVfaiqvrG7//hYagTWn84gwA6q6i+q6iVVdW1V/VlV3X98/suq6vyqemdVvauW3dIlAAADFElEQVSqHj8+/6+r6rVVdVmS/1lVdxp34t5TVW+qqsur6ruq6l9U1Rs2vc6jqur125Tw1CRv3KauE6rq7VX1L8fHPzau5bqqesn4aZckeUxVfcn4OScn+QdJ3jb+/u+Orw8MnDAIkNxtyzDxv9r0vY9090OTvDzJC8bnfjzJH3X31yf55iS/WFVfNv7eP0nyfd39yCRPTHJykjOSPHP8vST5oyQPqKoTx8fnJPmtber6xiTXbD5RVV+V5M1JfrK731xVj05yvyQPS3Jmkq+rqm/q7r9J8qdJzh7/6JOT/I++fdupq5M8YsL/f4B9zDAxwNGHiTc6dtdkFO6S5NFJHldVG+HwrknuO/76D7r7o+Ov/1mS13b3F5L8dVW9JUm6u6vq1UmeVlW/lVFIfPo2r/3VST686fjOSa5M8tzuvmpTLY9O8q7x8d0zCodvze1DxW8c/+8zNl3r1ow6hcDACYMAR/fZ8f/eltvfMyvJd3b3+zY/saq+IcmnNp86ynV/K8llST6TUWDc7v7CT2cUNDd8PqNQ+q1JNsJgJfm57v5v2/z87yb5pap6aJK7dfe1m7531/H1gYEzTAwwvSuS/HBVVZJU1UN2eN7bknzn+N7Br0pycOMb3X1LkluS/ESSC3b4+RuSnLrpuDPq7t1/04znK5I8o6ruPq7lpKq65/g1PpnkUJLzM+oSbvaPk1y/2y8K7H86gwDjewY3Hf9+dx9teZmfTvLLSa4bB8K/SPKYbZ73uiT/IqPQ9X+SvCPJkU3f/+9JTuzuP9/hdd6cUYD8w40T3X1bVT05yWVV9fHu/vWqekCSt4+z6SeTPC2jYeBkFAJfn9Ew8WbfPL4+MHB1+73EAMxaVd29uz9ZVX8/owkd39jdfz3+3q8leVd3v3KHn71bkreMf+a2Gdf11iSP7+6/neV1gfUjDALMUVUdSvIVSe6S5Be6+4Lx+Wsyur/wUd392aP8/LcmuaG7PzjDmk7MKGD+7qyuCawvYRAAYMBMIAEAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwP4fsYaj1JsgcyYAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 720x576 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "data091 = pd.read_csv(\"energy_2009_1.csv\")  \n",
    "x091 = data091['x']\n",
    "y091 = data091['y']\n",
    "xerr091 = data091['xerr']\n",
    "yerr091 = data091['yerr']\n",
    "\n",
    "\n",
    "plt.figure(figsize=(10,8))\n",
    "plt.loglog()\n",
    "plt.errorbar(x091, y091, yerr=yerr091, xerr=xerr091, fmt=',r')\n",
    "plt.xlabel(\"Energy (keV)\")\n",
    "plt.ylabel(r\"$keV^2 (Photons{\\ }cm^{-2}{\\ } s^{-1}{\\ } keV^{-1})$\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoMAAAHnCAYAAADD41dfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3X+0bvd8J/D3hyAp5qomwSRM0okhEYSqmhrtpaXRIVTLUMrE7ylaneoSq6Zd/TmlXaZLGzpURLRoEpSQNu2EG8OkSuKIaOhK0mrT6ELTXGIJEp/54zxHjuPce5/n3OfHec5+vdZ61n32fvaz9+fkx77v89n7+93V3QEAYJhutegCAABYHGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYsEMWXcAyOfzww/uYY45ZdBkAAAd0ySWXfKG7jzjQdsLgBI455ph89KMfXXQZAAAHVFWfGWc7l4kBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBAIABG3wYrKrjq+oPqurcqvpvi64HAGCe5hoGq+rQqvrrqvp4VX2yqn7lIPZ1RlV9rqou3+Szk6vq01V1ZVWdtr/9dPcV3f38JE9K8qCt1gMAsIzm3Rn8apJHdPf9k5yU5OSqesj6DarqyKq644Z1x22yrzOTnLxxZVXdOsnpSR6d5IQkT6mqE6rqvlX1ng2vI0ffOSXJB5NcePA/IgDA8phrGOxVN4wWbzN69YbNfjDJu6rq0CSpquckefUm+/pAkus2OcyDk1zZ3Vd399eSvC3J47r7E939mA2vz4329e7u/v4kT53GzwkAsCwOmfcBR527S5Icl+T07v7w+s+7+5yqOjbJ26rqnCTPTPLICQ5xVJJ/XLd8TZLv2089u5M8Icntkpy/j20em+Sxxx23WYMSAGB5zX0ASXff3N0nJTk6yYOr6sRNtnllkhuTvDbJKeu6ieOozQ67n3r2dPfPdPfzuvv0fWxzXnc/d9euXROUAQCw/S1sNHF3X59kTza/7+9hSU5M8s4kvzzhrq9Jcvd1y0cnuXZrVbLt7N69+gIApmLeo4mPqKo7jd4fluSHk3xqwzYPSPL6JI9LcmqSO1fVr09wmI8kuWdVHVtVt03y5CTvnkb9AAA7zbw7g3dL8v6quiyroe0vu/s9G7b5jiRP7O6ruvsbSZ6R5DMbd1RVb01ycZJ7VdU1VfWsJOnum5K8MMkFSa5IcnZ3f3JmPxEAwBKb6wCS7r4syQMOsM2HNix/Paudwo3bPWU/+zg/+xgMAgDALQb/BBIAgCGb+9QyMJGNg0VWVpK9e799/Z49cyoIAHYWnUGWjyl+AGBqdAbZ3jZ2/NY6gjqBADAVOoMAAAMmDAIADJgwCAAwYMIgAMCAGUDCcjFwBACmSmcQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDAhEEWb/fu1RcAMHfCIAAwGb/E7yjCIADAgAmDAAADJgwCAAzYIYsugIHZ7B6TlZVk797NP9uzZ8YFAcCw6QyyPezategKAGCQdAaZr806fWsdQV1AAJg7nUEAgAETBgEABsxlYgAYmoOdMHp/A//G5dagbUNnEACYnIF/O4bOIIvnt0OA+TrY866BfzuKziAAwIAJg+xMHqIOAGMRBgEABkwYBAAYMGEQAGDAhEEAgAEztQw7w8bBIptNiGoKBAD4NjqD7FwmRAWAA9IZZGfY2PUzISoAjEUYBAAm4xftHcVlYgCAARMGAQAGTBgEABgw9wyyM7mfBQDGojMIADBgwiAAwIAJgwAAAyYMAgAMmDAIADBgwiCzsXv3LY+EAwC2LWEQAGDAhEEAgAETBgEABkwYBAAYMI+j4+Dsa5DIykqyd++3f+4xcQCwregMMju7di26AgDgAHQGOTj76vStdQR1AgFgW9MZBAAYMGEQAGDAhEEAgAETBgEABswAEmbDwBEAWAo6gwAAAyYMAgAMmDAIADBgwiDDtXv3vh+nBwADIQwCAAyYMAgAMGDCIADAgAmDAAADJgwCAAyYJ5AwHBtHDq+sJHv3fuv6jU9OWfvME1UA2KF0Bhm2XbsWXQEALJTOIMOh6wcA30ZnEABgwIRBDsyTOgBgxxIGAQAGzD2DsN5WRhwDwBLTGYQDMeIYgB1MZ5Bb7Ou+wM26Y2uWuUu2We1GHAMwMDqDjEd3DAB2JJ1BbrGv7pfuGADsWDqDAAADJgwCAAyYMAgAMGDuGYT9cZ8kADucMMiBCUQAsGO5TAwAMGDCIADAgAmDAAADJgwCAAyYMAgAMGDCIADAgG0pDFbV7avq1tMuBgCA+RorDFbVrarqJ6vqvVX1uSSfSvLZqvpkVf12Vd1ztmXCNrV79+oLAJbUuJ3B9yf590leluSu3X337j4yycOS/FWS36qqp82oRgAAZmTcJ5D8cHd/fePK7r4uyduTvL2qbjPVygAAmLmxOoPrg2BVvfRA2wAAsBwO2BmsqrPXLyY5KckrZlYRAABzM85l4i9297PXFqrqtTOsBwCAORonDP7GhuVfnEUhsDTWjx5eWUn27v3WdXv2zLkgANi6A94z2N1/lyRVdfho+bpZFwVLZdeuRVcAAFs27mjiJDkjySmzKgSWxvrO31pHUDcQgCU1yRNIamZVAACwEJOEwZ5ZFQAALITOIADAgE0SBl82syoAAFiIscJgVd27uy+fdTEAAMzXuJ3B86vqjKq6x0yrAQAYx+7d3zrHK1s2bhi8d5KPJbmoqn63qo6YYU2wPPbsMa0MAEttrDDY3V/r7t9LcnySa5J8uKp+taruONPqAACYqUkGkKS7b+zu30ly3yQ3Jrm0ql4yk8oAAJi5icJgVR1TVScneXaSeyT5UpLfnEVhAADM3liPo6uqy5IcneQfknwqyRVJ3pfk9CSfnll1AADM1LjPJv6xJFd3t6eQAMAieSY6UzbuAJKrurur6p5V9Yaq+v1ZFwYAwOxNdM9gkjcnOTfJDyRJVZ1YVWdNvSoAAOZi0jB4q+7+syQ3J8noqSQnTr0qAADmYtIweG1VHZukk6SqKslhU68KAIC5GHcAyZoXJ/nDJHetqlOTnJzEM4sBAJbURGGwu/9+NM/g45PcP8lFSc6YRWEAAMzeRGGwqn69u1+e1UEk547W3XoWhQEAMHuTXiY+qqp+srvfkiRVdUSSs5M8fOqVAcBOszZH4MFYWUn27j24fZmjkHUmDYPPS3JBVV2Z1UEkb0zy0qlXBQDs265di66AHWTcx9GdleTSJB9L8oIkb0lyU5LHd/eVsysPAHaQaXTkdsoTSA62S6pDOjXjTi3zptG2z8xqEDwmyb8meVpV/cRsSgMA2A8d0qkYqzPY3RcmuXBtuaoOSXJCVkcUPySjwSTAOjvlt3eAWTjYc6Nz7NRMes9gkqS7b0py2ej15qlWBADA3Ez6BBIAAHYQYRAAYFK7d09nqqBtQBgEABiwicNgVZlXcFntoN9iAIDpOOAAkqo6e/1ikpOSvGJmFQEAMDfjjCb+Ync/e22hql47w3pg+a11XzebENUUCMDBch5hysa5TPwbG5Z/cRaFwI5kQlQAtrkDdga7+++SpKoO7+4vdPd1sy8Lltjab+0mRAVgCUwy6fQZSU6ZVSFMyf4GiBzoOY5CCwAMziSjiWtmVTA/LlsCAOtM0hnsmVXB9Oyvu+eyJQCwgc4gAMCATRIGXzazKgAAWIixLxN39+VVdVGSx3b3F6vq+UkOTfKa7v7azCoEANjILU9TM+nj6O40CoLfk+Q5Sb4zyeunXxbsAHv2OFkBsO1NMoAkSb5eVYckeXqSV3T32VX10RnUxSwIJgDABpOGwVcn+XhWLw+fNlp3h6lWBADA3EwUBrv7rKp6R5Kbu/srVXVckotnUxoAALM2aWcw3X3DuvdXJjl1qhUBADA3kw4gAQBgBxEGAQAGTBgEABiwA4bBqnpkVb2+qk4aLT939mUBADAP4wwg+emsDhJ5eVXdOclJsy0JAIB5Gecy8ee7+/rufkmSRyX53hnXBADAnIwTBt+79qa7T0ty1uzKAQBgng4YBrv7XUlSVYePln9v1kUBADAfk4wmPmNmVQAAsBCThMGaWRUAACzEJI+j65lVAQAwb7t3b/27KyvJ3r1b28eePVs/7gzoDAIAbMWuXYuuYCom6Qy+bGZVwBCt/Ta5zX5DBBiMgzn/7qBz+Nidwe6+PEmq6jazKwcAgHmapDOYqvrDJE+oqi8nuTbJZUkuM90MAMBymigMJnlYkrt099er6qgk909yv+mXNT9VdXySn01yeJILu/u1Cy4JAGBuJhlAkiR/leQ7k6S7/6m7z+/u3xr3y1V196p6f1VdUVWfrKqfnfD46/d1RlV9rqou3+Szk6vq01V1ZVWdtr/9dPcV3f38JE9K8qCt1gMAsIwmDYOvS3JRVb2kqh5WVZMOo7kpyc939/FJHpLkBVV1wvoNqurIqrrjhnXHbbKvM5OcvHFlVd06yelJHp3khCRPqaoTquq+VfWeDa8jR985JckHk1w44c8DALDUJg2Df5Tk7KxeXv7pJP+vqq4a98vd/dnuvnT0/ktJrkhy1IbNfjDJu6rq0CSpquckefUm+/pAkus2OcyDk1zZ3Vd399eSvC3J47r7E939mA2vz4329e7u/v4kTx33ZwEA2AkmvWfwmu7+5fUrqup2WzlwVR2T5AFJPrx+fXefU1XHJnlbVZ2T5JlJHjnBro9K8o/ra07yffupY3eSJyS5XZLz97HNY5M89rjjNmtQwgTWT066ccLSHTA9AQDLZ9LO4MrG+/y6+6uTHrSq7pDk7Ule3N1f3Ph5d78yyY1JXpvklO6+YZLdb7Jun09P6e493f0z3f287j59H9uc193P3bVDJpdkG/HfFAALNmln8C5JfriqXprk0iQfT7LS3eeMu4PRPIVvT/LH3f2OfWzzsCQnJnlnkl9O8sIJarwmyd3XLR+d1WlwYPHWd/90BAHYBibqDHb3k0aDP45N8ktJ/jb7uQS7UVVVkjckuaK7X7WPbR6Q5PVJHpfk1CR3rqpfn6DMjyS5Z1UdW1W3TfLkJO+e4PsAAIMxURisqgdV1e27+6vdfWl3v6m7XzLBLh6a5KeSPKKqVkavH92wzXckeWJ3X9Xd30jyjCSf2aSWtya5OMm9quqaqnpWknT3TVntJF6Q1QEqZ3f3Jyf5OQEAhmLSy8RnJXng2kJVHZ7kId39nnG+3N0fzOb39K3f5kMblr+e1U7hxu2esp99nJ99DAYBAOAWkw4gubG7b1xb6O4vJPnV6ZYEAMC8TBoGr66qR29Yd9tpFQMAwHxNepn4RUn+rKp+KquPprtPkrEnnQYAYHuZKAx292er6nuSPD7JSUkuS/LfZ1EY7HimlAFgG5i0M5juvjmr8wS+ffrlAAAwT5PeMwgAwA4iDAIADNikk05fVFX/ZvT++VX14tFTPgAAWEKTdgbv1N1fHA0ieU6S78wmE0IDALAcJh1A8vWqOiTJ05O8orvPrqqPzqAuAADmYNIw+OokH09yaJLTRuvuMNWKAACYm0nnGTyrqt6R5Obu/kpVHZfk4tmUBgDArG1lnsEb1r2/MsmpU60IAIC5mSgMVtUjkjw1yfVJLs/qE0gu7+6vzqA2AABmbNLO4B8lecHoe/fL6mPp7pPkuCnXBQDAHEwaBq/s7neO3p8z7WIAAJivSecZvKiqfq6qaibVAAAwV5N2Bu+T5MQkL62qS5KsJFnpbl1CAHa+3btX/9yzZ5FVsB3soP8GJp1a5glJUlWH5ZZg+JC4ZAwAsJQmHU185yQ/l+TIJH+T5KzuPnMGdQEAMAeT3jP4tiRfSnJeku9I8sGqevDUqwIAYC4mvWfwbt39ytH791TVnyR5S1YvFQMAsGQm7QxeV1X3W1vo7quz2iEEAGAJTdoZfF6Sc6vq/yb5RJITklw19aoAAJiLSTuDX0nywCTvz+ogko8necq0iwIAYD4m7Qy+s7sfmOTstRVV9ZAkfzXVqgBgVtbmCtyKlZVk796t72MHzU3HzjFWZ7CqnlRVv5XkjlV1fFXdet3Hr5tNaQCwDe3ategKYKrG7Qx+KMmhSZ6d5FVJ7lVV1ye5NquXjgFgORxMd84TSNiBxgqD3f1PSc6qqqu6+0PJNyegPjbJp2ZYHwAAMzTpPYNXVNWv5VufQPLl6ZcFAMA8eAIJAMCAeQIJAMCAeQIJAMtn9+6DmyIG+KaDfQLJfeIJJAAAS2uizmB3fyrf+gSSlXgCCQDA0pqoM1hVj0jy1CTXJ7k8yWVJegZ1AQAwB5NeJv6jJC8Yfe9+SR6f1UvFx025LgDYfkw2zQ40aRi8srvfOXp/zrSLAQBgvsZ9NvFZVfXiJBdX1c/PuCYAAOZk3AEkbxpte9ckP1VVn6mqd1fVr1XVE2dXHgAAszTus4kvTHLh2nJVHZLkhCT3T/LguGQMwFZsda7AlZVk796tfd99f/AtxgqDVfX0JK/KanfwvCQv7O7LsjqaGADmb9euRVcAO8K4A0h+Kckjk/xTkhcl+c3RnwCwdVvt0q11BHX54KCNGwa/2N0fG73/H1X14VkVBADA/IwbBu9WVc9NckWSTyW5zexKAgBgXsYNg7+c1Ummn5rkvknuUFXnJ/l4ksu6+60zqg8AgBkadzTx69YvV9XRWQ2H903yo0mEQQCAJTTuaOLq7m8+g7i7r0lyTZLz97UNAADb37iTTr+/ql5UVfdYv7KqbltVj6iqNyV5xvTLAwBglsa9Z/DkJM9M8taqOjbJ9UkOy2qY/Isk/6u7V2ZTIgAAszLuPYM3JnlNktdU1W2SHJ7kK919/SyLAwBgtsbtDH5Td389yWdnUAsAjMdk0zA1B7xnsKoeWVWvr6qTRsvPnX1ZAADMwzidwZ9OcmqSl1fVnZOcNNuSAACYl3FGE3++u6/v7pckeVSS751xTQAAzMk4YfC9a2+6+7QkZ82uHAAA5umAYbC731VV9163/HuzLQkAgHkZd9Lp86vqjI2TTgMAsNzGDYP3TvKxJBdV1e9W1REzrAkAgDkZKwx299dGl4ePz+oziT9cVb9aVXecaXUAAMzUuJ3BJKtPIunu30ly3yQ3Jrm0ql4yk8oAAJi5icJgVR1TVScneXaSeyT5UpLfnEVhAADM3liPo6uqy5IcneQfknwqyRVJ3pfk9CSfnll1AADM1LjPJv6xJFd3d8+yGAAA5mvcASRXdXdX1T2r6g1V9fuzLgwAgNmb6J7BJG9Ocm6SH0iSqjqxqjyRBABgSU0aBm/V3X+W5OYk6e7Lk5w49aoAAJiLScPgtVV1bJJOkqqqJIdNvSoAAOZi3AEka16c5A+T3LWqTk1ycpLLp14VAABzMWkYfHZWA+Djk9w/yUVJ3jTtogAAmI9Jw+BRSZ7Y3W9Ncm5VHZnkPUkePvXKAACYuUnD4POSXFBVVyX5RpIzk7x02kUBADAf4z6B5Kwklyb5WJIXJHlLkpuSPL67r5xdeQAAzNK4o4nfNNr2mVkNgsck+dckT6uqn5hNaQAAzNpYncHuvjDJhWvLVXVIkhOyOojkIVmdiBoAgCUz6T2DSZLuvinJZaPXm6daEQAAczPppNMs0u7dqy8AgCkRBgEABkwYBAAYMGEQgIPnNhZYWsIgAMCAbWk0MTO0v9+sV1aSvXv3v82ePVMuCADYyXQGl82uXYuuAADYQXQGt5v9dfbWOoK6f8BWOY8AG+gMAgAMmDAI241RmQDMkTAIADBgwiAAwIAZQLJM3PANzMNWblMYZ+qrzTivwcLpDAIwHaa+gqWkMwjAt9pKt86UNbC0dAYBAAZMZxBgp/FYS2ACwiBsF2t/OW/2l7W/nJkm9/YB6wiDsB35y5qD4bGWwASEQdgu1v5y9pc1AHNkAAkAwIAJgwAAA+YyMQAHz20NsLR0BgEABkwYBAAYMGEQAGDA3DMIMCTu7QM2EAZhu/GXNQBz5DIxAMCACYMAAAMmDAIADJgwCAAwYMIgAMCACYMAAAMmDAIADJgwCAAwYMIgAMCACYMAAAMmDAIADJgwCAAwYMIgAMCACYMAAAMmDAIADJgwCAAwYMIgAMCACYMAAAMmDAIADJgwCMtq9+7VFwAcBGEQAGDAhEEAgAETBgEABkwYBAAYMGEQAGDADll0AcAYNhs1vLKS7N27+Wd79sy4IAB2Cp1BWGa7di26AgCWnM4gLIPNOn1rHUFdQAAOgs4gAMCACYMAAAMmDAIADJgwCAAwYMIgwLzt3r35lEAAC2A0MSwro4gBmAKdQRgK3SgANiEMAgAMmDAIsEx0eIEpEwaB5SUYARw0YRAAYMCMJgaYlX11LVdWkr179/25keLAHOkMAizCrl2LrgAgic4g7EybdZz2143SiZqNff1zXft34J87sA3oDMKQ6EbNhoEswBLTGYSdaLOOk24UAJsQBgG2m/11GQ80+CQR+IGJCIPA8tgYgDYLRkMIQi73A1MkDALLbScGo/0FWpf7gSkTBoHlsTEAzTIYCV3AQBhNDAAwYDqDMBRD7nBtty7fdqkDIDqDAACDpjMIMIlJnu6iAwgsAZ1BgGnYiaOagUHQGQSYhKe7ADuMMAiwTAROYMqEQQD3AQIDJgwCy2t/l2zXG+d5vptxHyAwAMIgMAz7C3buAwQGTBgEdhbBDmAippYBABgwYRAAYMCEQQCAAXPPIMDBci8isMR0BgEABkwYBAAYMJeJgZ1vK5dxXfoFBkJnEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwIRBYH527159TWs7AA6aMAgAMGDCIPDtJu3M6eQBLC1hEABgwIRBAIABEwYBAAZMGAQAGDBhEABgwA5ZdAHADrZxhPHKSrJ377ev37NnTgUBsJEwCIwf2vZlku137dr68YVGgKkTBoHNbRbaJt1+Y3hbC3f7Wn8wxwdgS4RBYPzQti+Tbj/t4wOwZQaQAAAMmDAIADBgwiAAwIAJgwAAAyYMAgAMmDAIADBgppYBvt2spnQxVQzAtiMMAgdv2iFPaASYG5eJAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABgwYRAAYMCquxddw9Koqs8n+cyUd7sryd4p73MWFl3nPI8/q2NNe7/T2N/B7OPwJF84yOMzfYv+f3WeluVn3Q51zquGWR5nmvte9Pkzmc859N919xEH2kgYXLCqel13P3fRdRzIouuc5/Fndaxp73ca+zuYfVTVR7v7QQdzfKZv0f+vztOy/Kzboc551TDL40xz34s+f46+v23OoS4TL955iy5gTIuuc57Hn9Wxpr3faexv0f9emb4h/Ttdlp91O9Q5rxpmeZxp7tv5cx2dQWBLttNvtQDLZjudQ3UGga163aILAFhi2+YcqjMIADBgOoMAAAMmDAIADJgwCAAwYMIgMBVVdfuqelNVvb6qnrroegCWRVV9d1W9oarOXcTxhUFgn6rqjKr6XFVdvmH9yVX16aq6sqpOG61+QpJzu/s5SU6Ze7EA28gk58/uvrq7n7WYSoVBYP/OTHLy+hVVdeskpyd5dJITkjylqk5IcnSSfxxtdvMcawTYjs7M+OfPhRIGgX3q7g8kuW7D6gcnuXL0m+zXkrwtyeOSXJPVQJg4twADN+H5c6GcsIFJHZVbOoDJagg8Ksk7kvx4Vb02O+gxTQBTtOn5s6q+q6r+IMkDqupl8y7qkHkfEFh6tcm67u4vJzl13sUALJF9nT//Jcnz513MGp1BYFLXJLn7uuWjk1y7oFoAlsm2PH8Kg8CkPpLknlV1bFXdNsmTk7x7wTUBLINtef4UBoF9qqq3Jrk4yb2q6pqqelZ335TkhUkuSHJFkrO7+5OLrBNgu1mm82d196JrAABgQXQGAQAGTBgEABgwYRAAYMCEQQCAARMGAQAGTBgEABgwYRAYvKq6uapW1r1OW3RNa6rq3Kr67tH7Gyb43p6q+pEN615cVa+pqiOq6s+nXSuwnDybGCD5SnefNM0dVtUhowlmD2Yf90ly6+6+egtff2tWn25wwbp1T07yC939+ar6bFU9tLs/dDA1AstPZxBgH6rq76vqV6rq0qr6RFXde7T+9lV1RlV9pKo+VlWPG63/r1V1TlWdl+QvqupWo07cJ6vqPVV1flX9RFX9UFW9c91xHllV79ikhKcmedcmdR1eVRdX1X8eLf/CqJbLqupXRpudm+QxVXW70TbHJPm3ST44+vxPR/sHBk4YBEgO23CZ+L+s++wL3f3AJK9N8pLRul9M8r7u/t4kD0/y21V1+9Fn/zHJM7r7EUmekOSYJPdN8uzRZ0nyviTHV9URo+VTk7xxk7oemuSS9Suq6i5J3pvkl7r7vVX1qCT3TPLgJCcl+Z6q+oHu/pckf53k5NFXn5zkT/qWx059NMnDxvznA+xgLhMD7P8y8VrH7pKshrskeVSSU6pqLRwemuQeo/d/2d3Xjd7/pyTndPc3kvxzVb0/Sbq7q+rNSZ5WVW/Makh8+ibHvluSz69bvk2SC5O8oLsvWlfLo5J8bLR8h6yGww/klkvF7xr9+cx1+/pcVjuFwMAJgwD799XRnzfnlnNmJfnx7v70+g2r6vuSfHn9qv3s941JzktyY1YD42b3F34lq0FzzU1ZDaU/kmQtDFaS/9nd/3uT7/9pkldV1QOTHNbdl6777NDR/oGBc5kYYHIXJHlRVVWSVNUD9rHdB5P8+Ojewbsk2b32QXdfm+TaJC9PcuY+vn9FkuPWLXdWu3v3Xjfi+YIkz6yqO4xqOaqqjhwd44Yke5KckdUu4Xo2lAk+AAABDUlEQVT/IcnlB/pBgZ1PZxBgdM/guuU/7+79TS/za0l+N8llo0D490kes8l2b0/yQ1kNXX+b5MNJ9q77/I+THNHdf7OP47w3qwHy/6yt6O6bq+rJSc6rqi9292uq6vgkF4+y6Q1JnpbVy8DJagh8R1YvE6/38NH+gYGrW+4lBmDaquoO3X1DVX1XVgd0PLS7/3n02e8n+Vh3v2Ef3z0syftH37l5ynV9IMnjuvtfp7lfYPkIgwAzVFV7ktwpyW2TvLK7zxytvySr9xc+sru/up/v/0iSK7r7H6ZY0xFZDZh/Oq19AstLGAQAGDADSAAABkwYBAAYMGEQAGDAhEEAgAETBgEABkwYBAAYsP8PKU0Aa4bRWc8AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 720x576 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "data092 = pd.read_csv(\"energy_2009_2.csv\")  \n",
    "x092 = data092['x']\n",
    "y092 = data092['y']\n",
    "xerr092 = data092['xerr']\n",
    "yerr092 = data092['yerr']\n",
    "\n",
    "\n",
    "plt.figure(figsize=(10,8))\n",
    "plt.loglog()\n",
    "plt.errorbar(x092, y092, yerr=yerr092, xerr=xerr092, fmt=',r')\n",
    "plt.xlabel(\"Energy (keV)\")\n",
    "plt.ylabel(r\"$keV^2 (Photons{\\ }cm^{-2}{\\ } s^{-1}{\\ } keV^{-1})$\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnIAAAHnCAYAAAA8bbD4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzs3Xt81OWZ///XnQQkWCSCYCOHooVqVMBDWrXUKmJApSoLrGCt8Iu2Wndb7dKq+LBV3H6rsa1drWCrVRGwFZQqrAW0tlTXsgINhoIFFBUKGYOEQxCWYA7cvz8mM5mZfGbm88mcJ+/n4+HDzCczkzuBMNdc931dl7HWIiIiIiK5pyDTCxARERGRzlEgJyIiIpKjFMiJiIiI5CgFciIiIiI5SoGciIiISI5SICciIiKSoxTIiYiIiOQoBXIiIiIiOUqBnIiIiEiOKsr0AtLlhBNOsEOGDMn0MkRERETiWrdu3R5rbb949+sygdyQIUOorq7O9DJERERE4jLG/NPN/bS1KiIiIpKjFMiJiIiI5CgFciIiIiI5qsuckRMREZH0aG5upra2liNHjmR6KVmvR48eDBw4kG7dunXq8QrkREREJKlqa2vp1asXQ4YMwRiT6eVkLWste/fupba2lpNPPrlTz6GtVREREUmqI0eO0LdvXwVxcRhj6Nu3b0KZSwVyIiIiknQK4txJ9OekQE5EREQybsrjbzHl8bcyvYyco0BORERE8s7OnTsZPXo0ZWVlnHHGGTzyyCMA7Nu3j4qKCoYNG0ZFRQX79+8H/OfVbr31VoYOHcqIESN4++23g881b948hg0bxrBhw5g3b15Gvp9oFMiJiIhIRi2p8VGzo4E12/YxqmolS2p8CT9nUVERDz30EJs3b2b16tXMmTOHTZs2UVVVxZgxY9i6dStjxoyhqqoKgBUrVrB161a2bt3KE088wS233AL4A7/77ruPNWvWsHbtWu67775g8JcNFMiJiIhIxiyp8XHXixtpaj0KgK+hkbte3JhwMFdaWso555wDQK9evSgrK8Pn87F06VKmT58OwPTp01myZAkAS5cuZdq0aRhjOP/882loaKCuro5XX32ViooK+vTpw/HHH09FRQWvvPJKQmtLJrUfERERkbQLnIer2dEQDOICGptbuWPxBp5bu4NFN1+Q8Nfavn07NTU1nHfeeXz88ceUlpYC/mBv9+7dAPh8PgYNGhR8zMCBA/H5fFGvZwtl5ERERCRjIoO4eNe9OnToEJMmTeLhhx/muOOOi3o/a22Ha8aYqNezhQI5ERERSbtFN1/AopsvYEBJsePnB5QUJ5yNa25uZtKkSVx33XVMnDgRgBNPPJG6ujoA6urq6N+/P+DPtO3cuTP42NraWk466aSo17OFAjkRERHJmNvHnUpxt8Kwa8XdCrl93KkJPa+1lhtvvJGysjJmzJgRvH7VVVcFK0/nzZvH1VdfHbw+f/58rLWsXr2a3r17U1payrhx4/jjH//I/v372b9/P3/84x8ZN25cQmtLJp2RExERkYyZcPYAAO5YvIGm1qMMKCnm9nGnBq931qpVq1iwYAHDhw/nrLPOAuD+++9n5syZXHPNNTz11FMMHjyYF154AYArrriC5cuXM3ToUHr27MncuXMB6NOnDz/60Y/44he/CMA999xDnz59ElpbMhmnvd98VF5ebqurqzO9jC5vSY2Pn736Lh81NHJSkn5ZRUQku2zevJmysjJPjwkUPySjuCHXOP28jDHrrLXl8R6rjJykTaDEvLG5FWgvMQcUzImIdHFdMYBLBgVyknLpLDEXERHpSlTsIGmT6hJzERGRrkYZOUlYvHNvgUzbqKqV+BoaOzw+GSXmIiIiXZEycpKQwLk3X0MjltijVVJVYi4iItJVKSMnneL23Bu0Z+QCWTpVrYqISAdzx/v/X7kss+vIMcrISUK8nnubcPYAVs28hG1V41k18xIFcSIikhI7d+5k9OjRlJWVccYZZ/DII48AsG/fPioqKhg2bBgVFRXs378f8DcQvvXWWxk6dCgjRozg7bffDj7Xjh07GDt2LGVlZZx++uls3749E9+SIwVy4mhJjY9RVSs5eeYyRlWt7LBV6na0is6+iYiIKy1NsGsDHPw4KU9XVFTEQw89xObNm1m9ejVz5sxh06ZNVFVVMWbMGLZu3cqYMWOoqqoCYMWKFWzdupWtW7fyxBNPcMsttwSfa9q0adx+++1s3ryZtWvXBsd6ZQMFctKBzr2JiEjaHdgBn34CbzyYlKcrLS3lnHPOAaBXr16UlZXh8/lYunQp06dPB2D69OksWbIEgKVLlzJt2jSMMZx//vk0NDRQV1fHpk2baGlpoaKiAoDPfOYz9OzZMylrTAadkROg/cwbJPnc28FdsLgSJj8DvU5M/TciIiK5IXAmDvzZuEO7/B+vmwu7NkJhd//tJJyZ2759OzU1NZx33nl8/PHHlJaWAv5gb/fu3QD4fD4GDRoUfMzAgQPx+XzU1tZSUlLCxIkT2bZtG5deeilVVVUUFhY6fq10UyAnHXTm3FvUs25v/BR2rPa/w/raL5K1RBERyScHdrR/bC007IC+Q5Py1IcOHWLSpEk8/PDDHHfccVHv5zSy1BhDS0sLb775JjU1NQwePJgpU6bwzDPPcOONNyZlfYlSICdA+GiUhPq9Rb7D8q31fxz5DgtUmSQi0pUFXgMO7oJHRoZ8wsKRBpj8dMI7Oc3NzUyaNInrrruOiRMnAnDiiSdSV1dHaWkpdXV1wfNuAwcOZOfOncHH1tbWctJJJ9Hc3MzZZ5/NKaecAsCECRNYvXp11gRyOiMnHcQ993ZwF8y9PP6B1AM7AOP/OPAOS0REJNQbPwUbseNjjyZ8Vs5ay4033khZWRkzZswIXr/qqquYN28eAPPmzePqq68OXp8/fz7WWlavXk3v3r0pLS3li1/8Ivv376e+vh6AlStXcvrppye0tmRSRk46iHvuLdZ2aYd3WIFUtYt3WDpPJyLS9dSuhdam8GutTf7rCVi1ahULFixg+PDhnHXWWQDcf//9zJw5k2uuuYannnqKwYMH88ILLwBwxRVXsHz5coYOHUrPnj2ZO3cuAIWFhfz85z9nzJgxWGs599xz+da3vpXQ2pLJOO0J56Py8nJbXV2d6WVkRLwRWnFF2y41BTCg3Hm79A8zoGZB+C9nYXc4+/roZ+X+MMO/BXtupc7TiYjksM2bN1NWVubtQV24IbDTz8sYs85aWx7vscrI5blAK5HG5lagvZUI0LlmvMHtUhv7QKqbd1huztN1wV9oEZEuSf/ed4oCuTwT2kYE3LUSAWIXMXRmu/Tbf/W2cLcBooiIiASp2CHPeW0lElOyD6RWLvP/N/kpf0DoFCCKiIhIVMrI5ZnIzFpCrUQipehAaswAUWflREREolIgl+duH3dq2Bk5cBih5bZa1Ot2qVtRAsT31q1k3F+Xda5AQ0REpAtQIJfn4rYSgcxPX4gIEJNeoCEiIlmv8pVKAOZeNjfDK8ktCuS6gLARWnPHw3r8/4G/WvSjav9WptP0BUhbJVGgUCNWgYYCORERcWPnzp1MmzaNXbt2UVBQwE033cRtt93Gvn37mDJlCtu3b2fIkCE8//zzHH/88Vhrue2221i+fDk9e/bkmWee4ZxzzgHgzjvvZNky/2vhj370I6ZMmRL8OpMnT+anP/0pp5xyCnfffTfz589n//79HDp0KHif2bNnc+yxx1JZWZn071PFDl3dgR3+KlHImukLSS3QEBGRrHbxoosZPm841R9XU/1xNcPnDWf4vOFcvOjihJ63qKiIhx56iM2bN7N69WrmzJnDpk2bqKqqYsyYMWzdupUxY8ZQVVUFwIoVK9i6dStbt27liSee4JZbbgFg2bJlvP3226xfv541a9bws5/9jE8++QSAf/zjH7S2tgbHd1155ZWsXdvx3PgNN9zAL3/5y4S+n6jfZ0qeVbJXaHatM9MXUihQfBGrQCNUwo2ORUQk4/Ye2evpululpaWUlpYC0KtXL8rKyvD5fCxdupTXX38dgOnTp3PxxRfz4IMPsnTpUqZNm4YxhvPPP5+Ghgbq6urYtGkTF110EUVFRRQVFTFy5EheeeUVrrnmGn77298GR3wBnH/++Y5r6dmzJ0OGDGHt2rV86UtfSuj7iqSMXFeWovl2iYo765X2c3S+hkYs7efoltT40rxaERHpjMpXKoPn4mLdJxm2b99OTU0N5513Hh9//HEwwCstLWX37t0A+Hw+Bg0aFHzMwIED8fl8jBw5khUrVnD48GH27NnDX/7yF3bu3An4x4Cde+65rtZQXl7Om2++mZTvJ5QycvnGy7zSVLUTSVCsAg2doxMRES8OHTrEpEmTePjhhznuuOOi3s9pZKkxhrFjx/K3v/2NL3/5y/Tr148LLriAoiJ/+FRXV0e/fv1craN///5s2bKlc99EDArk8o2XCtRUtRNJgrACDQcxz9F5CWZFRCQjAtWpw+cNj3ufzmpubmbSpElcd911TJw4EYATTzyRuro6SktLqauro3///oA/AxfItAHU1tZy0kknAXD33Xdz9913A/D1r3+dYcOGAVBcXMyRI0dcreXIkSMUFxfHv6NHCuRyVeic0oBo80pD5fgsu0VfPxkWV3Jl72+y8UCPDp8fUFIcP5hVoCcikvestdx4442UlZUxY8aM4PWrrrqKefPmMXPmTObNmxc843bVVVcxe/Zspk6dypo1a+jduzelpaW0trbS0NBA37592bBhAxs2bGDs2LEAlJWV8f777zNkyJC463nvvfcYNWpU0r9PnZHLJ8F5pWRNBWrStQVpj5T+seM5OtPMD7s/B9VPgT1Ka/UzXPnDX3PyzD8wqmpl+/m50EBPREQyqm+Pvp6uu7Vq1SoWLFjAypUrOeusszjrrLNYvnw5M2fO5LXXXmPYsGG89tprzJw5E4ArrriCU045haFDh/Ktb32Lxx57DPBn9S688EJOP/10brrpJp599tng1ur48eODhRMAd9xxBwMHDuTw4cMMHDiQWbNmha3n0ksvTeh7cmKc9oTzUXl5ua2urs70MlInUIHaEpLiLeoBt23I7axTaOYxNONoCljRazL/b98lfNTSi5OKDnL7CauZ0PIKHPoYsHxqi1jUejH3tNwA+AO9X/RbxuUHF/uLOsKe47j2s3hDC5WxExFJwObNmykrK/P0mFxsCNzY2Mjo0aNZtWoVhYWFUe9XU1PDL37xCxYsWOD4eaeflzFmnbW2PN4atLWaL7rCvNJgxtGCtVx+9H+4/PMftW8XHxwNjywk0E7lGNPCvxb+D4+2TKSeEhptNw7t24UtOooBWi3s37cHX4v/8Gug8nXE51/ilB2r+fD393D9rilqbyIikga5FMAFFBcXc9999+Hz+Rg8eHDU++3Zs4cf//jHKVmDtlbzRZZWoCascpn/v8lP+XvcOfW8C3AIZgs4yneLXgSgH/u5svCtwOYzhRzlXwr/Sj8agvf/TPMeBmz7PdijnLTt9zQ11Km9iYiIRDVu3LiYQRxARUWFq3N0naGMXL7I4grUpHCTcXQIZo8xLZxbsBWAW4teosCEHyUIBHqB7ddbi14iECwabNjn1N5ERESyjQI5yQ1uMo4hwWygYXBjc2vwWnnh+3SnJewpQgO9fuznXwvf4BjTEvxc6NYsQO/WvTD3cp2fExGRrKBATnKDx4yjU1Phd8ctpywkmxYZ7N1a9BKG2Bm7u4592X2fPhERkRRTIJdDNFvUm3hNhSODvfO6fcAxNnrGbrCp5+rWPwJHw/v05XhvPhGRbPDP66cB8LkF8zO8ktyiYoccodmiqTHh7AGsmnkJ26rG84V718OsA8H/lly9iVE9XuJrTQ8woOgTnun1GIVtGbtWC7/fUczJ794S3qNORESyws6dOxk9ejRlZWWcccYZPPLIIwDs27ePiooKhg0bRkVFBfv37wf8DYRvvfVWhg4dyogRI3j77beDz3XnnXdy5plncuaZZ7Jo0aKwrzN58mQ+/PBDDh8+zPjx4znttNM444wzgv3pAGbPns3cuampylVGLksFZooGxJot+tza9sa/i26+IC3r6wrCMnoHd8Ej3yNQCFHIUa4wb1HFFHwNhrte3Bh8jIiIeGebmvj0gw9oqa+nyOX80liKiop46KGHOOecczh48CDnnnsuFRUVPPPMM4wZM4aZM2dSVVVFVVUVDz74ICtWrGDr1q1s3bqVNWvWcMstt7BmzRqWLVvG22+/zfr16/n000+56KKLuPzyyznuuOP4xz/+QWtrK6eccgqHDx/mBz/4AaNHj6apqYkxY8awYsUKLr/8cm644QZGjRpFZWVlEn5S4ZSRyxExZ4tK6sVpbdLY3MrPFr/hL4Q4+HEmVigiktOaPvqIo4cOUT/nsaQ8X2lpKeeccw4AvXr1oqysDJ/Px9KlS5k+fToA06dPZ8mSJQAsXbqUadOmYYzh/PPPp6Ghgbq6OjZt2sRFF11EUVERxx57LCNHjuSVV14B4Le//W1wxFfPnj0ZPXo0AN27d+ecc86htrY2+LkhQ4awdm3yW4IpI5elIjNro6pW4mto7HC/ASXFysKlQ5zWJgDfNr9XIYSIiEuBM3GHq6v9YyXbNCxcSMPChWAMPcvLk3Jmbvv27dTU1HDeeefx8ccfU1paCviDvd27dwPg8/kYNGhQ8DEDBw7E5/MxcuRI7rvvPmbMmMHhw4f5y1/+wumnnw74x25de+21Hb5eQ0MDL7/8MrfddlvwWnl5OW+++SZf+tKXEv5+QimQyxG3jzu1QzuN4m6F3D7u1AyuqgsJqZp1Cqr7sZ9rCt/wZ+2cCiEO7nIe+xXtuohIF1E8YgRNO3fSum+f/0JBAYUlJXSP02TXrUOHDjFp0iQefvhhjjvuuKj3cxpZaoxh7Nix/O1vf+PLX/4y/fr144ILLgjOWq2rq6NfxDZwS0sL1157LbfeeiunnHJK8Hr//v3ZsmVLUr6nUNpazVYHd4Vt0004ewAPTBzOgJJiDP5M3AMTh+tMVgbcPu5UiruFz9SbUfR7itr6z2EtNOwIf9AbP23P1rm5LiKS5z63YD6fWzCfIYsW0mtshf+iMWAtvcaOZcjC5xLOxjU3NzNp0iSuu+46Jk6cCMCJJ55IXV0d4A/E+vfvD/gzcDt37gw+tra2lpNOOgmAu+++m/Xr1/Paa69hrWXYsGGAf0TXkSNHQr8kN910E8OGDeN73/te2PUjR45QXFyc0PfjRBm5bBX6At+2TRevnYakR2TbkuG9G7mmeRWFR9vHhx05tI8L93yd7j98jh/2Wcnlnzzv/1Rotq6lCXxrw69/8zX/bWXqRKQLadmzl8J+/ejWvz89Roygpb4+4ee01nLjjTdSVlbGjBkzgtevuuoq5s2bx8yZM5k3b17wjNtVV13F7NmzmTp1KmvWrKF3796UlpbS2tpKQ0MDffv2ZcOGDWzYsIGxY8cCUFZWxvvvvx8cv/XDH/6QAwcO8OSTT3ZYz3vvvceoUaMS/r4iKZDLpLnjna87vcAXdg+/j3qXZVRYUP2HGVATnpI3tr2RcMO+PbQWFVDIUVotLNlRzA+abuKh7o9zdYH/eocsnkMgLyKSrwbNfjR4Zq703nuS8pyrVq1iwYIFDB8+nLPOOguA+++/n5kzZ3LNNdfw1FNPMXjwYF544QUArrjiCpYvX87QoUPp2bNnsF1Ic3MzF154IQDHHXcczz77bHBrdfz48bz++utceuml1NbW8pOf/ITTTjstWGTxne98h29+85vB9dx7771J+d5CKZDLRgd2AAaw7S/wfYdmelUSTYxCiH7s518K3/QHa7S3LfkNl3GFeSt4HSwcaYAnK/wdTqJl6kRExJWvfOUrjufeAP785z93uGaMYc6cOR2u9+jRg02bNjk+z+TJkxk9ejT33XcfAwcOjPr1ampqOOOMMzjhhBM8fAfuKJDLJKes2sFd8MhIoH2bjiMNMPlpbbFlq5BCiCEzw/9Mf1z0tOPYr4e7z+5wHXu049k6p/N2IiJ5KBcnOhQXF3Pffffh8/kYHKM4Y8+ePfz4xz9OyRpU7JBtHPqVYY/qMHyOGFASfpD1nIKtHGM6jv36nNnd4TqtTVBc4g/cg9oCefWmExHJSuPGjYsZxAFUVFQEz9ElmwK5bOOwTUdrk/+6ZL3IitbxTQ9Q1rooOO5ryJHfMeTI7yj7dF7w41E9XmofDfa5UbED+YhqZhGRbBVtm1HCJfpz0tZqtgnZppPcE1nRelJJMbePOzV4PW4vwHiBvIogRCQH9OjRg71799K3b1+MMZleTtay1rJ371569OjR6ecwXSViLi8vt9XV1ZlehnRxS2p8jkFerOt3LN5Ac2srw4tqeanoLn+BhCmAAeVqVyIiWam5uZna2toOPdakox49ejBw4EC6desWdt0Ys85aWx7v8crIiaSRUy/AJTW+sEydr6GRu17cSPU/9/H7db62ebqGf+U1WmwBhUbtSkQku3Xr1o2TTz4508voEpSRE8mwL9y9oi1Yi64f+3nzmO/RwzS3XyzqAZ8d4S9w/qjaf5bOFMCMLcrKiYjkOLcZORU7iGRYvCAO4Nail6K3Kzmwo71AwloVRoiIdCEK5EQyLLJlSUBhyAFhpzYmUduVrP+tP3jTHFcRkbynQE4kwyJbloC/mvXa8wYFr49vesDftqStlUmsdiWfNrew7GfTaK2e6//curn+iREiIpJ3VOwgkmGxWpaUf65Ph+sAo6pW8lFDI68Wr+QLtuN4sPN5hxZrKDSEF0aoulVEJK+o2EEkh0RWuII/e9fcepSWo/7f5aiFEbdt8G+zrpsL51aGV7cqwBMRySoqdhDJQz979d2wIA6gsbk1GMRBjMKIX43yB3FO2606TycikpO0tSqSQz5qaIx7n6iFEY37w6tbG3bA3PHQ0gS+tskR6+bCro3tjYZFRCSrKZATySEnlRTjcwjmSoq78WnLURqbWxnf9ADg33KddO4A/rKlnuaGj3izx39wTPAR1l/t2trkb18SvBzRaFhERLKatlZFcki0CtdZV53BAxOHM6CkGIO/pcmkcwfw+3U+fA2NfLfopQ7Vrdij0OfzHduXHGlQ7zkRkRyhjJxIDolV4Rr6efBXtgbO00Xbbt3795fpZVroHnrdHm0f9aUiCBGRrJaTgZwxZgIwHugPzLHW/jHDS3It2nB0Ebec5rU6CT1PF9huDSjuVkhjcyvLut9FX/PP8Ae2NkFt25k5zXAVEclqaW8/Yox5GvgasNtae2bI9cuAR4BC4ElrbZWL5zoe+Lm19sZ4982G9iPRWkc8MHG4gjlJulFVKx3P00UzoKSYVTMv8d+ILILQDFcRkbRy234kExm5Z4DZwPzABWNMITAHqABqgb8ZY/4bf1D3QMTjb7DW7m77+Idtj8sqUx5/y/F6zY6GDnM1G5tbuWPxBp5bG37AfNHNF6RsfdI13D7uVMc3DpHtSwI6VMQe2AEYwLbPcNV2q4hIVkl7sYO19n+AfRGXvwS8b6390FrbBCwErrbWbrTWfi3iv93G70FghbX27XR/D50VbTi6m6HpIl5NOHtAhwKIwG0nJ4Ven/xUWxFEIGOvGa4iItkoW87IDQB2htyuBc6Lcf/vApcCvY0xQ621v3a6kzHmJuAmgMGDBydpqfFFy6ZF2+oaUFKsDJykRLTzdE6ZutGn9eMLd6+gqfUovzh2PlcfbSWsPrblU39T4cP1/tvr5sJFdyorJyKSQdnSfsQ4XIt6eM9a+0tr7bnW2m9HC+La7veEtbbcWlver1+/pCw0EdFaRwTmZ4qkg1OmLtCqJJAdPrV5M4W2OeKR1t9UOHjTtmflDu6CuZerbYmISJplS0auFhgUcnsg8FGG1pIy8VpHiKRLZKYukIkLCFS5hhVAHNwFj4yElkAbk7bt1ovuVHWriEiGZEsg9zdgmDHmZMAHTAW+ntklpYbb1hEi6RTtnGZYAcQbP+3YVFjbrSIiGZX2rVVjzHPAW8CpxphaY8yN1toW4DvAq8Bm4Hlr7T/SvTaRrspVAUTtWn+PuTDabhURyaS095HLlGzoIyeSraL1OJx07gCe/1stTa1HGeB0FCC43Xqk/VpRD7htgz+gWzcXzq3UdquIiEfZ3EdORLKM0/nN0af1CyuA8DU0cteLG4OPuWPxBn5knmRKUcSIr8B2a+Ne/1ZsYLsVq/5zIiJJpkBORICO5zdDZ7UGNDa3Muu//8GnLUdpaj3KOd230p2IGa6B7dbAebrQ7VYVRIiIJFW2tB8RkSzTYdJDm4bG5mCAN77pAYYc+R1DjvyOUT1eglkH4PvvQmHoe0QL1U9D9VPtGbonK/yfijxHp3N1IiKeKJATEUcnRSmAiCYY+DlVt4a2hbQWGna03zd0SoTT1AgFdyIiUSmQExFH0RpYH9+zm+P9g4GfY3VrKOsf//X4ReFZusjbgcBNI8FERKJS1aqIRLWkxtehgTU4j/iKWuH6hxlQsyA8uCvsDqYQWgLbt8Zf7Rp6+zMnQu/B4FvbdqkAZmxRoYSIdAmqWhWRhMVqYO22wnWCU4autYnwyXw2JIhru/1/u8O3aANFEyqUEBEJUkZORBIWOeIrIGzEVyinLF2kgm7+QM6GVM4W9YAb/wSv3KE2JiKS19xm5PL+jJwx5kpjzBMHDhzI9FJE8parEV/4t2q/cPcK/rH2T3HO0QFHm8ODOPAHdi9+M/zMnIohRKQLy/utVWvty8DL5eXl38r0WkTy1YCSYnwO7UpCK18D0yOaWo8yvvUBwH+27oGJw523b3/9Fdi1MfxaaxPUb/F/HGg0HFoMoW1XEeli8j4jJyKpF63CNVAcAf4zdU4Nhn/26rvOT/rtv/r70oX+95nPtn/eWv8EichKV2XoRKQLyfuMnIikntOIr8i5rNEaDEe73sHBXf62JUEWDteH3NQECRHpehTIiUhSxKpwBf82a7zt15gcGw2HapsgEWg+HNh6VUGEiOQxba2KSFq42X6NKW6jYegwQeKNB7XVKtJZ+t3JCQrkRCQlAhWqQ2YuY1TVSgAemDicASXFGPwFElELHZxEnpn77PA4D7Cw/rfwp1nRJ0M9QQEtAAAgAElEQVTohUokuscvgn/+r6aqZDn1kRORpAtUqEZOf/AUuHkVrzedKYCb/ie8B90fZvi3YM+t1Hk6kYC546GlqeNUFSwsrlQPxzRRHzkRyRjPFarJEG/r1drwHnS/qVDFq0g0B3a0fxw4pqAMXVZSRk5Eku7kmctw+pfFANuqxqdnEQd3wSMjoeWIw0IKoOcJ/jFggZWV3+D/UBk66eocf3cMwTOomnucFsrIiUjGRKtEdV2h6mBJjY9RVSs5ue3M3ZIaX+wHxKpytUfh8N7QC/6K18gMnRvK4km+cfzdcSgkkqygQE5Ekq4zFaqxArXAmTtfQyMW8DU0cteLG2MHc3G3WlsjL4R86OGFKnSyhEg+iFsh3lZIpDcvWUFbqyKSEktqfDEbBEfeN1pxxHNrd1Czo8Fxnmv3wgLe+8nl8RcTrxDCSeExUDoCpvzWv4V0cFf7Qe/Aoe/LfgqPf8V/f203Sb5y+v0p7A5nXx/9CELo74t+JzrF7daqGgKLSErEaxAMMOXxtwAcA7XG5lbuWLyBsweXOAZxQNTrHUTLMHx2uL+tidMLVWsz1P6tfUJEZOZtx2p/8URAIIt30R0dAz6nFzO90EmucPr9aW3yX4/m8Yvg0C5NWEkDBXIiknGxArVFN1/AqKqVjlMhBrg9c/ftv8b+vGOg17amdXPBtw7q1vtvh06PqN8Scv+27abmwx0DPqcXs9DAUC90ks3i/f6ECrQuObTLf1sTVlJOZ+REJGMW3XwBi26+IGpAFrie8FSIeCKbDZff2P45a6E+tG1KjOMoLUfg78/5D4pHFk/UbWwvightfaKzRpJvnFqXSMrkfSBnjLnSGPPEgQMHMr0UEYkiXqA24ewBcadCeK5qjebgLn9wFWShpWM2ML6I4omnKtp7cIW90B0Nf6FLtApWVbSSSZOfgiMNIRdUGJFqeR/IWWtfttbe1Lt370wvRUSicBOoTTh7AKtmXsK2qvGsmnlJhyDOc1VrNLHaljgp7O7P4H3/XSjqEeVOIcFg9dPQuL/9U61N4S90sapg3QRpqqKVdPpNBTwwKPzvb+TvT+SbFUkqnZETkazgpjjCyZTH34pZLOH5OeO2XogQOPTtOgC0cLQ5/NLRVnj8q3DdYv+WKzifLYp2ri5QONF0uP0s3/rf6mySpN6BHfDpJ+1/JztTGCEJUSAnIjkv0arWJTU+7li8gabWowwo+U9uvzp6q5Sofv0V9wFgZMB3tNl/ODxWFWysIC1QIVgUctYwkAVRIYWkglNRw66N3gojJCnyfmtVRPKbm2KJWALbsoGgr9PbspEFE6GFE4Xdw+8b2I6ddQAGfKn9ulMV7K8v9J+tCy24CN2q+k1F+4tp6Fm+yC1bkWSLLGpo2BH9vpIyCuREJOclUtV6x+INYY2IoX1bNinibTUdiPHi13KkfR5stCAt1uN1NklSxamo4UiD3jhkgLZWRSTnBbZB3U6SCJVws+F4Ym01HdwV8WLogT0Kr90b+/E6myRudKY5dayiBm3np5UCORHJC50tlhhQUpxYs+FEOL0YBkYfXXQHPDLSn5Vz0toEG5+HggLnx0d7Mf1NBezZAt+pbn/RjvdCrikU+a0zUxhU1JA1tLUqIl3CkhofX7h7BUMi+sylvNlwLLFeDKMFeYGzdZ/5LNhW/ygxp8dDe7uS0GbEoVWGAY9f1N7jzqnFSejnJX/MHR9+xnLdXHiywt1jo50JVbFD2ikjJyJ5L1pBAyS2LZuwWC96TlWwgSAt9MW3qAfctsE5UxbItDxV4T9j99q97Y9b/1v/6LG970PTIf+1dXNh88v+c3mB7MzBXeGPUUuT/KKChZxnrI0xbiaPlJeX2+rq6kwvQ0TSLFqfOYDuhQW895PLM7CqBP381PbgKtpW6m8qwBdjm6uwO5jC6FMrAgHiGw+297aLt20r2S90mxzbcfs+1hsDSStjzDprbXm8+ykjJyJ5L5GChvAec87ZOjf3SZrIAolABWtkpixWNWvgcbEEiik2vRT/a0nuCD0PBypYyAM6IycieS2RPnNueswlrQ+dW25GICVSDRsQKKaw4a1Z1NIki8Ub4RZ5Hu6f/6uChTygQE5E8l5nCxpi9Zib8vhbTHn8rdT3oYvkplrQ67zYaOIVU0TjFFC4mRPbGal63lwUrygl8jzc576cnQULXv5MI2e9dkHaWhWRvNfZggY3W7Ip70MXyc2LbLx5scV9/AUObkeKld/YvtUWaF9y8OP27dXAtT6nwL4P/eesQgsmoHMtLtzw+rxO7VeSJfK5U/m1QjmNy4rc/u6QpbXZu03u5c80ctZrF6RATkS6hM70mYvVY27RzRcAMKpqZeb60EUTL9j79Vf8czEjfXZ4+2NDiyVCX/CdXjgD1+rf9RdPfPpJ+OOwHatl930I/99yeOWOzvevi9wqdBOUpPKFP/K50xlkRGbbQquOF1dCyeey9zzc3PHtHzvNb40ccef2vpXL4n/tdAXbKaStVRGRKNxsyWa0D11nuekBFhYYtL3ghwZOgRFh0Wa9Bh732r3w+FfDr9W/6w9wnqqIvRXodavQ6X6h23SRgV8yt+NC27Ssm9ueVUrF14rkNC4r8OcT+Bm+92punIfz0g4lGa1TnPoq5hi1HxERiWFJjS/ulqyb++SUg7uc21Icc1z77FeAgm5QfHz4tc5wankRmhE0BTBjS8etQqc13vin8CxfoFVL+Y2wZVl7cIWB8huSl436w4z2Ni0Y/1qCgW2Ur5WsbNAfZkDNgohAzUDPE+BwfdtNh59hJsT6nqP9mTq1Q3Fz39BMn5OWpvC/YwPKo2f/wF2GL4ncth9RRk5EJIYJZw9g1cxL2FY1nlUzL3EM0NzcJ6c4FUscbYXDeyOuNXe81hlOlbDRsm2BDNufZjlvFYZm+UIzcNVPQ+P+0Du3Z60SdXCX/7lCnzssOxnla7nNBsUrHnE8E2nDv99oGct0i/U9u6nI7sx9Y60l+NjcbYac92fkjDFXAlcOHTo000sREckNToHB0Wbn+0a2J+mMQH+6cyv92bTLfhr9YH5gy3L3FuetwoB1c/0ZqdDniPwe7FH41Sho/TR6Vixe1uzgLv/WcbyfQ+h5tFjFCU5fz+nwf+g1pzORwYxVS/v3n6nihkBmLNqZtkCmy8v8Vjf3jZVBC/x8gqz/79zkp6P/fALfR5ozc/Foa1VEROKLViARkwFCXmNMIRQUdGxpAuGTJvqd5i+G8LpVGDrxIvD14gVYpghsS3hlbqjQrdlYn3ejqBgKu/m/v4Ydztu8oV+v/l3n7T9L7G1ncN5uzdRkjkAAtPf98O/5MydC36GZCYw68/NJcyCnyQ4iIpI8gayP23NMTi+UthVaowRWofer3+JwhyhbhYEXXacmyJFBXOQLdeg5PKeqV6eqWGx7Je3Cb7R/3s0ZtECQ1ufzsCu0z6D1P79vXfjXG1Aef/sv8ucQ4CW7lWqVy2JnwDIhm34+CVIgJyIi7sU6mxQaTETrZRfa4iRg7viO2ZrI4oB4W4VumiBHvlA7ncOr29i+ten0+UDBxBsPRm/54SQ0KNz4vD8zGcpafwYu9LZTwBd2zs/h5xCQ6aa+kdz+vUmXbPv5JECBnIhImniZyZrW+a1euM1keHmhnPxUx2xNZHASLxDwEjhC9Aa5xxznP4z/2r0dP1/9NMGt4uqnIyoc45xBCwv6nDKTNjzLiXUO+Fqb/TvWYQ/Nkn5wseRRBizbKJATEUmDaDNZAcd2Jm7vm3apyGS4ydbECwS8rsvpa7YcaQ+mNjznP88Wvqjwj52KJ5wCKqdt38jtaNdb0UfDlwEdAqJsfBMwxfwMSgk20pbkUSAnIpIiUx5/K/hxzY6GDmO7AjNZn1sbfu7JzX1jvSAGvm7OvGi6ydYkO4CMN8YMnIsyQkUGgtEyTJ0NVCF6RjGKrH4TICmhQE5EJA28zGRNZH7rkhpfMBAcVbUyK7IxcWXivFLk14xXxJFIFaiHQDU0+A9yuhaFlzcMbiTjzUBO/p100tLUcc5wFuhUIGeMORY4Ym0yGgiJiOSn0BfBWDNZI18svdw3VLKyMdm4NZdynTmD5/aMVxoD1UTeBKRCXmUI0zk71wNXfeSMMQXAVOA64IvAp8AxQD2wHHjCWrs1hetMmPrIiUgmBV7QGpvb3/8WdyvkgYnDo56Ri3Vfp8yNUzYGoHthAWcPLgnedhMMullnXonWJ8/j1mamxXoTsGrmJZ6fzzFD6IHbv5NuZPSowMFd8FDb/ORoY8OSLNl95P4C/Am4C3jHWv/bFmNMH2A0UGWMecla+2xnFywiks8CQZCbmaxe7hvKbTYm1otzZ7bmcuYsXiw5FKzFcvu4Ux0D8dvHnZqR9WRbhrDT3vhp+8dZViXsNiPXzVob89Snm/tkkjJyIpLv3GZjYgVya7bti/q5807u43g9LwK5PJJNW+PJzhBmhNsm2EmW1IxcaIBmjLnTWtthKm02B3EiIl2B22xMrMCrs+fzJHtMOHtA1myDZ1uGsFOyrZlxhLiBnDHm+dCbwFlAh0BOREQyq7NbsqHy4oVXskbg7162ZAg7JcubGcfdWjXGPGmt/WbI7V9Za29J+cqSTFurIiLuZNPWnEjWmDve///KZWn5csncWv1JxO27O7ckERHJBdm0NScisRXEu4O1dhuAMeaEttvRT8KKiIiISNrEDeRCPJ2yVYiIiIiIZ14COZOyVYiIiIiIZ15GdMVvOCciIpIFVLAhXYWXQE4ZORERyXqZmu+p4FEywdVkBwBjzJnW2ndSvJ6kM8ZcCVw5dOjQb23dmtXjYEVEuoTAZAkvDYa9zPxMZL5nZ5sed9kZtZIybtuPuB3RdZq1dktSVpYh6iMnIpJ5nc1aeQnkOjNmLFHJGg6v6RkSkNQRXcByY8zrwCxrrfPUZBERkRgS2fL0EuAkMmbMS8AYKm+Gw0vOcZuR6w7cDMwAlgI/sdbWp3htSaWMnIhIargNftK15ZmJbc68GA4vWcVtRs5V+xFrbZO19lGgDKgF1hhj/tMY0yvBdYqISBeRrqzVhLMH8MDE4XQv9L/EDSgpTvlZtdvHnUpxt8Kwa5pRK+ngutgh7EHGHAvcBlQCj1trf57shSWbMnIiIpmV71krVa1KMiX7jFzgSYcApwGnAoOBg8D9QNYHciIiklm3jzvVccszX7JWmlErmeAqkDPGbAAGAjuALcBmYCUwB3g3ZasTEZG8EQhyfvbqu3zU0MhJylqJJMxtRu5fgA9tZ/ZhRURE2ihrJZJcbosdPrDWWmPMMGPMU8aY2alemIiIiIjE5iqQC7EAWAx8FfzTHowx85O+KhERERGJy2sgV2CtXQG0ArSN7Doz6asSERERkbi8BnIfGWNOBiyAMcYAxUlflYiIiIjE5an9CPA94Engs8aYSuAy4J2kr0pERERE4vIUyFlrtxtjLgMmACOBN4CnU7EwEREREYnNa0Pg/2et/SH+gofFbdcKYz9KRERERFLB6xm5AcaYrwduGGP6AX9K7pJERERExA2vZ+RuBl41xryPv+BhLnBn0lclIiIiInG5HdE1H3gbqAH+Hfgd0AJMsNa+n7rliYiIiEg0brdW57Xd9wb8QdwQYD/wDWPM5NQsTURERERicZWRs9b+Gfhz4LYxpgg4HX/l6vm0FT6IiIiISPp4PSMHgLW2BdjQ9t+CpK5IRERERFzxWrUqIiIiIllCgZyIiIhIjlIgJyIiIpKjPAdyxpic6htnjLnSGPPEgQMHMr0UERERkaQy1trYdzDm+dCbwFnW2mEpXVUKlJeX2+rq6kwvQ0RERCQuY8w6a215vPu5qVr9xFr7zZAn/lVCKxMRERGRpHCztfqTiNt3p2IhIiIiIuJN3EDOWrsNwBhzQtvtfalelIiIiIjE56XY4emUrUJEREQkxaY8/hZTHn8r08tIKi+BnEnZKkRERETEMy+BXOzyVhERERFJK2XkRERERHKUl0DurpStQkREREQ8cx3IWWvfMca8YYw5DsAY821jzPeMMd1TtzwRERERicbriK4Sa+0nxphzgW8BxwO/Sf6yRERERJJnSY2Pmh0NrNm2j1FVK1lS48v0kpLCzWSHUM3GmCJgGvCgtfZ5Y4zmXomIiEjWWlLj464XN9LUehQAX0Mjd724EYAJZw/I5NIS5jWQ+yXwd6AHMLPt2meSuiIRERGRCIn0f6vZ0RAM4gIam1u5Y/EGnlu7w/XzLLr5gk6vIVU8BXLW2vnGmBeBVmttozFmKJBfnfVEREQkr0QGcfGu5xKvGTmstYdCPn4fqEzqikREREQiJJING1W1El9DY4frA0qKszLL5oXXYgcRERGRnHL7uFMp7lYYdq24WyG3jzs1QytKHs8ZOREREZFcEihouGPxBppajzKgpJjbx52a84UOoEBOREREuoAJZw8IFjbk+nZqqLhbq8aYCmPMb4wxZ7Xdvin1yxIRERGReNxk5P4Nf0HDD40xfYCzUrskEREREXHDTbFDvbW2wVr7A2As8MUUr0lEREREXHATyC0LfGCtnQnMT91yRERERMStuIGctXYpgDHmhLbbj6Z6USIiIiISn5eq1aeBq1K1EBEREZFUyqdq1QAvDYFNylYhIiIiIp55CeRsylYhIiIiIp4pIyciIiKSo7wEcnelbBUiIiIi4pnrQM5a+w6AMaZb6pYjIiIiIm55ychhjHkS+NgYs9MYs6ZtdNd3U7S2pDDGXGmMeeLAgQOZXoqIiIhIUnkK5IALgROttYOAicBLwLFJX1USWWtfttbe1Lt370wvRURERCSpvPSRA1gNHA/sttb6AB+wPOmrEhEREZG4vGbkngDeMMb8wBhzoTFGaS4RERGRDPEayD0LPI8/k/dvwP8aYz5I+qpEREREJC6vW6u11tp7Qy8YY45J4npERERExCWvGbn1xpjbQi9Yaz9N4npERERExCWvGbkTgUuNMXcCbwN/B9Zba19I+spEREREJCZPgZy19hoIbqeeAQwHzgMUyImIiIikmadAzhhTDmy21v4f/ozc2ylZlYiIiIjE5fWM3HygNXDDGHOCMeZryV2SiIiIiLjhNZA7Yq09Erhhrd0D/GdylyQiIiIibngN5D40xlweca17shYjIiIiIu55rVr9LrDCGHM9/nFdZwBqCCwiIiKSAZ4yctbaOuBc4PdAP2AD8PUUrEtERERE4vCakcNa24o/kPt98pcjIiIiIm55PSMnIiIiIllCgZyIiIiIg8pXKql8pTLTy4jJUyBnjHnDGHNc28ffNsZ8zxijqlURERGRDPCakSux1n5ijDkX+BZwPPCb5C9LREREROLxWuzQbIwpAqYBD1prnzfGVKdgXSIiIiISh9dA7pfA34EewMy2a59J6opERERExBVPgZy1dr4x5kWg1VrbaIwZCryVmqWJiIiISCyd6SN3KOTj94HsLucQERERyVOeAjljzCXAdUAD8A7+yQ7vWGs/TcHaRERERCQGrxm5Z4F/b3vcCGAC/nmrQ5O8LhERERGJw2sg97619qW2j19I9mJEREREMu3iRRez98je4O3h84YD0LdHX16f8nqGVuXMax+5N4wx/2GMMSlZjYiIiEiGhQZxbq5nkteM3BnAmcCdxph1wHpgvbVW2TkRERHJadk+jsuJ1/YjEwGMMcW0B3Xno21WERERkbTzWrXaB/gPoD+wCZhvrX0mBesSERERSau5l80F2s/E5QKvZ+QWAgeBl4GewF+NMV9K+qpEREREJC6vZ+RKrbU/bfv4D8aYRcDv8G+vioiIiOS8vj36OhY29O3RNwOric1rILfPGDPCWrsBwFr7oTGmZwrWJSIiIpIRgRYjgeKHwJZrNvIayN0MLDbGvAlsBE4HPkj6qkREREQkLq9n5BqBc4C/4C94+DtwbbIXJSIiIiLxec3IvWStPQd4PnDBGHM+sDqpqxIRERGRuFxl5Iwx1xhjqoBexpgyY0xhyKefSM3SRERERCQWtxm5VUAP4JvAL4BTjTENwEf4t1tFREREJM1cZeSstT5r7Xzgamvt5dbaU4BLgXuBS1K5wEQZY640xjxx4MCBTC9F2jTv3s32b1xPS319ppciIiKS07wWO2w2xvzYGPM4cD3wobX2/1KwrqSx1r5srb2pd+/emV6KtNk2aTKN1dXUz3ks00sRERHJaV6LHRYCfwLWAMPxT3aotNauTfrKJG/88/ppAByurgZrg9cbFi6kYeFCMIae5eV8bsH8TC1RREQkJ2myg6RN8YgRNO3cSeu+ff4LBQUUlpTQffDgzC5MRKSLuXjRxVEnFwSa4Up2NwIO0GQHSbnQTFvdrFk0LHoe0707tqmJXmPHUjrr3gyuTkSk63EK4mJdl+yV6GSHM9BkB/GgZc9eSqZO5fgp17B/0fMqeBARSaLASKlUPUcuZKi6Gk+BnLV2izHmHGAC/vFc64Hvp2Jhkp8GzX40+HHpvfdkcCUiIiK5z1MgZ4y5BLgOaADeATYANuaDREREJC3cZsyGzxue8HNIdvC6tfos8O9tjxuBPzN3BjA0yesSERERkTi8BnLvW2tfavv4hWQvRkRERFKvb4++UatWJbe4CuSMMfOBt4G3jDHft9Y+lNpliYiISKqoxUj+cDvZYV7bfT8LXG+M+acx5r/bpjz8a+qWJyIiIiLRuMrIWWv/DPw5cNsYU4S/anUk8CW0zSoiIiKSdq4ycsaYacaYPcaYfcaYeUCxtXaDtXaBtfb2FK9RckDz7t1s/8b16gsnIiKSRm63Vu8BKoDTgB3A/SlbkeSkbZMm01hdTf2cxzK9FBERiWH7lKm8W/5FvfHOE26rVj+x1ta0ffwjY8yaVC1IcsM/r58GwOHqarDtrQQbFi6kYeFCMIae5eVA+IguERHJrKaPPuLooUPUz3lMIxLzgNtArtQYcxOwGdgCdEvdkiSXFI8YQdPOnbTu2+e/UFBAYUkJ3QcPTvi5m3fvxjfj+wz8r19Q1K9fws8nIpIPAm+kvXLzxtstvUHPHm4DuXvxNwC+DhgOfMYYsxz4O7DBWvtcitYnWSr0l7hu1iwaFj2P6d4d29REr7Fjk/Iub89jv6Jx3Tq9axQRSYJUvvGWzHFbtfpE6G1jzED8gd1w4ApAgVwX1rJnLyVTp3L8lGvYv+j5Tp+7cLNdW7Z5U1LWLCKSqxLJhtXNmkXDwkVgDFibtDfekjluGwIba9tfWa21tUAtsDzafaTrGDT70eDHpffek/Dz6V2jiEhqtOzZS2G/fnTr358eI0ao4CEPGDexlzHmdeD3wFJr7Y6Q692BrwDTgb9Ya59JzTITV15ebqurqzO9DHEp7F0jUDJlit41iohIl2GMWWetjXt40e0ZucuAG4DnjDEnAw1AMf72JX8E/stau76zixWJ1LJnLyXXXpvwdq2IiGTOxYsujjrTVWPCksNVRi7sAcZ0A04AGq21DSlZVQooI+edqkZFRCQRw+cNj/q5jdM3pnEluSfZGbkga20zUNepVUlO2TZpMq319aoaFRERACpfqUz7c829bG7SvmY+ihvIGWMqgGuAOdba9caYmyKrWCX3hfYlUpNfERGR3OAmI/dvQCXwQ2NMH+Cs1C5JMk1VoyIi4sRrdizW1qoybcnhJpCrbzsL9wNjTBXwxRSvSTIgMsOWqia/yaZzfCIi0pUVuLjPssAH1tqZgPbUuoBAk98hixZSMnUqLXv2ZHpJjrZNmkxjdTX1cx7L9FJERCRC3x59PV0X79z2kTvNWrslDetJGVWt5odo0x+CIqY/KGMnIiK5yG3VqpuMHMByY8zTxhgdkspRzbt3s/0b1+dNP7biESMo7NOn/UJBAYV9+lA8cmTY/ZSxExGRfOY2kDsNqAHeMMY8bIxRaiPH5EtA87kF8/ncgvkMWbSQXmMr/BdDZgYOWfgc/7x+GpvLTmfzaWW0tgWuDQsXsvm0MjaXnR5WoSsiIpLLXPWRs9Y2AY8aY34DfAdYY4x5FviZtfZgKhco3kQGKW5aiUButhOJNf1BlbciItIVeJ7sAGCMORa4DX9bksettT9P9sKSrauckYsM5GxTU9SAxnTrFrxfLgZy8Wheq4iI5KqUTHYwxgzBv816KjAYOAjcD2R9INdVOAVkudJKJNk0r1VERPKdq0DOGLMBGAjsALYAm4GVwBzg3ZStTpIi0EqkqwU0g2Y/Gvy49N57MrgSERGR1HDbfuTzwIe2M/uwWaKrbK2Kd2pRIiIi2SapW6vW2g/annQYMBNotNZ+J7ElimSHbZMm01pfT/2cx7rElrOIZL+LF13M3iN7O1zv26Mvr095Pf0Lkqzltv1IwAJgMfBVAGPMmcaY/DslL3kvVouSLSM1TlhEMsspiIt1XcJVvlJJ5SuVmV5GWngqdgAKrLUrjDH3A1hr3zHGnJmCdYmkXLQWJacsXZLZhYlI3kpGcBHrOTSIvuvxmpH7yBhzMmABjDEGKE76qkRSLFZTYZ2TExGRXOE1I/c94Engs8aYSuAy4J2kr0okTdSiRETSyW3GbPi84Qk/RzYL9DzNxx6m6eY1kPsm/uBtAjASeAOYl+xFiaSLWpSISDptnzKVTz/4gM+vWN6ls/+2qYlPP/iAlvr6pP4cIotEAgFxPheJeN1aHQD8q7V2sbX2R/gLH/6Q/GWJiIjkn6aPPuLooUNx51737dHX0/Vc4/bn4FVXLBLxNKLLGNMdeBW4EzgKPAPcaa1dlpLVJZH6yImISCIiRyB6ETn3Oihi7rUb2bAd2dmfRWd+Dlv2belwbeF3z3C8b/XH0V/ny0+M/3POpm1rt33kXGXkjDHzjTHfAy4A/h14Avg1MCEXgjiRzmrevZvt37heZ+dEJCHFI0ZQ2KdP+4WCAgr79KF45MjMLSoD9HNIPreTHcbgPxMX+G8I8DdgFfCOtXZxCteYFPmSkdMUgvR678Kv0lpfT8nUqWoWLDqOLvcAABPtSURBVCJhvDbtrZs1i4aFi/xV8kDJlCld8t+VVP4cYhWJbJy+MfhxLhRbJHuyw5+BP4c8eRFwOv6g7nz8Z+UkDTSFIPX+ef20Dun/hoULaVi4EIyhbPOmDK5ORLKF1/NYqpL3a9mzl8J+/ejWvz89Rozosj+HZPFatQqAtbYF2ND234KkrkiAjucPYgUWoecKsvndRS6J1iy4++DBmV2YiKRVZxv4Oj7uazD3Mn+lfFeukh80+9Hga1yyfw59e/SNmiXNV50K5CT9FFikTyAYDkv/tzULVhZURCR7Bba0A4F0NhUvpIoCuSzllFmrmzWLhkXPY7p3xzY1KbBIMW2DiEisQCDfm/amknaPkkeBXA5p2bOXkqlTFVikiZoFi4hItlMgl0MUWIiIZI+ueB4rV8TLiKZqskQmKJCTlPNaoi8ikgrJ/rdI/37lrtDJErl+REmBnKRcvo5MUU8/kdySr/8WdUXJmiwRrQNELNl2vk+BnKSMm7L9ylcqc/ZQsHr6iWSfpLYLQUUL+SYfO0AokBPxaHPZ6WoWLCKSQYlkxfKttVTeB3LGmCuBK4cOHZrppXQ5gXey+Vain4/v6ETyhdqFSDz5NlmiINMLSDVr7cvW2pt69+6d6aVInhiyaCG9xlb4b4S8oxuy8LnMLkxEROIaNPtRjhkyhIKePSm9956wjhC5KO8DOcm8aKX4uVyiH2gWfPJLL1IydSote/ZkekkiEkc+/lskYmzIWZ98Vl5ebqurqzO9DBEREcmwQNVrtlWghjLGrLPWxi2lVUZOREREJEflfbGDpI4a/YqIiGSWAjnpNDXXFJFcs33KVD794AM+v2K5Gnl3Ydm8peqVAjnxxG2zzc40+lWGT0RSLZ9GM4mAAjnJIsrwiUioaG/uSj4t4tk3o/eEc5KPo5lEQIGceBSaZUtWc818H+UlIp0T7U1cwzEtnp9LjbwlXymQExGRtGrevRvfjO8z8L9+wbfW3dGp55h1XaHj9Vhv+PJtNJMIKJCTBPTt0TfqmTYv8nWUl4g42zZpMq319dTPeQzOT9/XzbfRTCKghsCSRWIFchunb0zjSkTEjUBTVbciz6kFRTmn9rVLa6I+l/5NkHynhsCSczQ+RyS/FY8YQWGfPu0XCgoo7NOH4pEjM7cokRynrVXJGmoxIpJbnKo447UR8nJOrW+M5xIRPwVyEqQ+biKSqHhthLycU9O/OyLxKZCTIPVxExG33DYH7/CYrwFf+3zblW0qZBJJkAK5Ls7LpAbIjepRZRZFRKSrUCAneUeZRZHUi/amTm2ERNJLgVwXl4pJDakWK+N2cu+TM7AiERGRzFD7Eck5yriJZC+1ERJJL2XkJChZkxpSxc15vrmXzY2ZWRSR1NI5VJH0UiCXAaFzBov69cv0coK6yj/AKoYQEZF8oa3VDNg2aTKN1dX+OYPi2tzL5ro6qxdva0dbsyIdNe/ezfZvXK/5oyI5Rhm5FIg2fzByzmDDwoU0LFzoOGfQqWO6F1056xTr++tM7yuRriB0kH20SQsikn0UyKVR8YgRNO3cSeu+ff4LBQUUlpTQffDgpH+tfM46ZftZPoDtU6by6Qcf8PkVy7Nq+1xyX7y/W4kOso/1BjNSom84RSRxCuRSINY/bnWzZtGw6HlM9+7YpqaYcwa98JJpCr1vNrUVcSuRjGK6iiGaPvqIo4cOKbshSZfsv1vpfIMpIsmnQC7NWvbspWTqVI6fcg37Fz2v8yh5ZnPZ6Y7ZDXPMMZz29/UZXJnkEqesmtvMWWeyZF4G2YtIdlEgl2aDZj8a/Lj03nuS9ryRmbVcae6bCancmo2W3Thl6ZKEn1u6tlRmzrwMsheR7KJATrqcVBZ7DFm00DG7oXNy4kW0rFqqMmehbzBFJLcokMtTuVAQkK0SrfhVdkNSpWXPXkquvVZHM0QkyNiQMxf5rLy83FZXV2d6GZIDYm1Lb5y+MY0rERGRrsoYs85aG7t0HDUEFgmjPnMiIpJLtLUqIiJA124kLpKrlJETCdHVK3qla8vnRuIi+UoZOZEMat69G9+M7zPwv36hylZJi84eH4h8nN70iGQHZeREIkSr7E1Fxe+2SZNprK6mfs5jSX9uERHJf6paFUmDyE79kV36g0K69GuOpaSbKrZFsoeqVkWyWPGIERT26dN+oaCAwj59KB45MnOLEhGRnKMzciJp4JRd60yXfp2pyw7Z+OeQjIpTNRIXyT0K5EQypDNd+rdNmkxrfT31cx7TUPMMSuTPIVVBYDIqTtViRCT36IycSJYKPVcX70xdrp6n2z5lKp9+8AGfX7E8azJboTpzthFin29878Kv0lpfT8nUqQkF45FVpNUfR//3rfzE9rWp2lQkN7g9I6eMnEgOKB4xgqadO2ndt89/oaCAwpISug8enNmFJajpo484euhQxjOMbrNkbv8cIgNA6BgENixcSMPChR2CQFChi4i4p4ycSI4IO1MHlEyZknPbq4EAJ1Zmq2zzpjSvyluWzM2fg1MgZ5uaogaBplu3sPt2JpBTxalIflFGTiTPtOzZS2G/fnTr358eI0a4OlOXrTqTYfS6DRvt/rG2rCOzZE4BlZs/h2iBWGcKXEREYlEgJ5IjBs1+NNNLSFhogOM1qPG6Devm/p0JKBP5c0hlMK6KU5GuSVurWU5DrCVf7fzOdynq1y+sanfQ7Ec7XWDg9v6R2bJ82LIWkfyjrdU8oSHWkq9CM1ul994T9X5es2Ze759PW9Yi0vUokMsyXgZah95XLQUkXySjebKX++fDlrWIdF0K5EQk63ltntyZZssiIrlIZ+SynFoKSGelaoKA2+fNxjFWuUJnY0XE7Rm5gnQsRkTSb9ukyTRWV1M/57G4990+ZSrvln/R9ZgwN8/r5etLOJ2NFRG3lJHLcnpnLm4kOs4rVkNcN018Q6tHs63Zby7ozLgtnYsVyW+qWs0TCtbEq3hVm14a4np5Xq/3ExGRxCkjJ5KHYvVGCw3k3IyNitrEl+g919SbLZzXzLrOxoqIMnIiKZArW92xeqPFbIgbp1WH255r6s0WTmfeRCRVFMiJeJArL8heeqN5adXh9nnVm83PTV/IylcqO5x307gtEXFLgZxIF+d2woKkTzZld0UkuymQExFJkUCmLdaZN1Wfikgi1EdOREREJEcpIyci0klui1905k1EUkWBnIgHekGWUG6LX3TmTURSRYGciAd6QRbofDWqiEiy6YyciIiISI5SRk5ExCNVo4pItlBGTkRERCRHKZATEemkaEUuKn4RkXTR1qqISCep+EVEMk0ZOREREZEcpYyciEgEt41+RUQyTYGciHQJXoIzt41+RUQyTVurItIluAnOKl+pjNvs100zYBGRdFFGTkTymtspDCIiuUiBnMj/3979xlpW1Wcc/z5FBQLJmMBoKohowIKNCdoR02Ab1ADTlD8JGB2ViKBOTKrRF5pgNDbGFyRt0vhCRccwM9q02IqITCGOCo6EhigDGAEJhRBbRiTgn2AwVuP464uzr3N6c++dOzP7nH3WOd9PQuaedfZe+3eZzL7PXWuvdaSOG/1Kao1BTtJcM5xJmmc+IydJy7jRr6RWOCInaSGccMwJq65aXc4tRiS1wiAnaSEYziTNI4OcpOa5ga+kReUzcpKa5wa+khaVQU5S09wDTtIiM8hJkiQ1ymfkpJ74nNYwdmzeseYecZI0zwxyUk98TmtyDMmStDKnViXNvIOFZDfwlbSoHJGTNNPWs5jBUTlJi8ogJ2lQTptK0uEzyE2BP6ik1R1s2tTFDJK0Op+RmwIfgl8MPqd16A42beoecZK0NkfkenQ4P3SWn7Nj846+ytGUObo6OYfygfeStEgMcpIGc7Bp06VfbAzJkrSyJoNckjOBDwAnArdV1bUDlwSsPpq2nh9UkiRJh2rqz8gl2Z7kqSQPLGvfnOThJI8muXqtPqrqoap6L/BmYNMk65U0WT5bKEmHb4gRuZ3Ap4EvLTUkOQr4DHAesA+4O8nNwFHANcvOv6qqnkpyMXB119dM8/kerWfl8qKubp7n702SJm3qQa6q7khy6rLms4FHq+oxgCRfBi6pqmuAC1fp52bg5iS3AP+60jFJtgJbAU455ZRe6j8c/qDSelYut7a62XAqScOblWfkTgIeH3u9D3jtagcnORe4FDgauHW146pqG7ANYNOmTdVHoVLf1rPaeemYWXqmch7DqSS1ZlaCXFZoWzV4VdUeYM+kipG0svVusbPecDpLwVSSWjQrQW4f8OKx1ycDTwxUizRVS2HmUFY3O2UpSYLZCXJ3A6cneSnwE2AL8LZhS5Jm11BTluOBcj3B0613JGmyph7kklwPnAucmGQf8PdVdV2S9wG7Ga1U3V5VD067NmlS1rNyeT3HtPo8nSRpMoZYtfrWVdpvZY2FC1LL1jPd2dqUaF/hVJJ0+FK1GIs5N23aVHv37h26DKkXa01Z3n/F/VOsRJI0CUnuqaqDfujBrDwjJ2kGuIhCktoy9Y/oknTkJvWxVu77JkltcUROalDfo2Pu+yZJbTLISXPO6VJJml8GOWnOrWe61H3fJKlNBjlpTrnnnCTNPxc7SPqjSS2ikCRNxtzvI5fkIuCi00477T2PPPLI0OVIU+eec5LUnvXuIzf3I3JVtauqtm7YsGHoUiRJkno190FOWnROl0rS/HKxgzTn3GJEkuaXI3KSJEmNMshJkiQ1yiAnSZLUKIOcJElSowxykiRJjTLISZIkNcogJ0mS1CiDnCRJUqMMcpIkSY0yyEmSJDXKICdJktQog5wkSVKj5j7IJbkoybZnnnlm6FIkSZJ6NfdBrqp2VdXWDRs2DF2KJElSr+Y+yEmSJM0rg5wkSVKjDHKSJEmNSlUNXcNUJHka+O8Jdb8BaGE1xdB1TvP6k7pW3/320d+R9HEi8LMjvL76NfS/02lq4XudhRqnVcMkr9Nn34ty33xJVW082EELE+QmKcm2qto6dB0HM3Sd07z+pK7Vd7999HckfSTZW1WbjuT66tfQ/06nqYXvdRZqnFYNk7xOn3173/z/nFrtx66hC1inoeuc5vUnda2+++2jv6H/XtWvRfr7bOF7nYUap1XDJK/TZ9/eN8c4IictsFn7zVKSZt2s3TcdkZMW27ahC5CkxszUfdMROUmSpEY5IidJktQog5wkSVKjDHKSJEmNMshJ+qMkxyX5YpIvJHn70PVI0qxL8rIk1yW5YYjrG+SkOZdke5KnkjywrH1zkoeTPJrk6q75UuCGqnoPcPHUi5WkGXAo982qeqyq3jVMpQY5aRHsBDaPNyQ5CvgM8DfAK4C3JnkFcDLweHfY/inWKEmzZCfrv28OyiAnzbmqugP4xbLms4FHu98kfwd8GbgE2McozIH3B0kL6hDvm4PyRi0tppM4MPIGowB3EnAjcFmSa5mjj7CRpB6seN9MckKSzwGvSvKRaRf1nGlfUNJMyAptVVW/Bq6cdjGS1IDV7ps/B9477WKWOCInLaZ9wIvHXp8MPDFQLZLUgpm8bxrkpMV0N3B6kpcmeR6wBbh54JokaZbN5H3TICfNuSTXA3cBf5ZkX5J3VdXvgfcBu4GHgH+vqgeHrFOSZkVL981U1dA1SJIk6TA4IidJktQog5wkSVKjDHKSJEmNMshJkiQ1yiAnSZLUKIOcJElSowxykpqWZH+SH4z9d/XQNS1JckOSl3VfP3sI5+1JcsGytg8m+WySjUm+0XetktrkZ61Kat1vquqsPjtM8pxu888j6ePPgaOq6rHDOP16RrvG7x5r2wJ8uKqeTvLTJOdU1X8eSY2S2ueInKS5lOTHST6R5N4k9yc5o2s/Lsn2JHcnuS/JJV37O5N8Jcku4JtJ/qQbAXswyX8kuTXJm5K8McnXxq5zXpIbVyjh7cDXV6jrxCR3Jfnb7vWHu1p+mOQT3WE3ABcmObo75lTgRcCd3fs3df1LWnAGOUmtO3bZ1Opbxt77WVW9GrgW+FDX9lHg9qp6DfB64B+THNe995fAFVX1BuBS4FTglcC7u/cAbgfOTLKxe30lsGOFus4B7hlvSPJC4Bbg41V1S5LzgdOBs4GzgL9I8tdV9XPg+8Dm7tQtwL/VgY/i2Qv81Tr//0iaY06tSmrdWlOrSyNl9zAKZgDnAxcnWQp2xwCndF9/q6p+0X39OuArVfUH4Mkk3wGoqkryz8DlSXYwCnjvWOHafwo8Pfb6ucBtwN9V1XfHajkfuK97fTyjYHcHB6ZXv979edVYX08xGqGTtOAMcpLm2W+7P/dz4H4X4LKqenj8wCSvBX493rRGvzuAXcD/Mgp7Kz1P9xtGIXHJ7xkFyguApSAX4Jqq+vwK598E/FOSVwPHVtW9Y+8d0/UvacE5tSpp0ewG3p8kAEletcpxdwKXdc/KvRA4d+mNqnoCeAL4GLBzlfMfAk4be12MRtXOGFtZuxu4KsnxXS0nJXlBd41ngT3Adkajc+NeDjxwsG9U0vxzRE5S645N8oOx19+oqrW2IPkk8Cngh12Y+zFw4QrHfRV4I6PA9F/A94Bnxt7/F2BjVf1olevcwij8fXupoar2J9kC7Eryq6r6bJIzgbu6XPkscDmjqVMYBbgbGU2tjnt917+kBZcDz85KksYlOb6qnk1yAqPFB+dU1ZPde58G7quq61Y591jgO905+3uu6w7gkqr6ZZ/9SmqPQU6SVpFkD/B84HnAP1TVzq79HkbP051XVb9d4/wLgIeq6n96rGkjo3B4U199SmqXQU6SJKlRLnaQJElqlEFOkiSpUQY5SZKkRhnkJEmSGmWQkyRJapRBTpIkqVH/Bxgt7r7jIkWRAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 720x576 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(10,8))\n",
    "plt.loglog()\n",
    "plt.errorbar(x00, y00, yerr=yerr00, xerr=xerr00, fmt='o', label='2000')\n",
    "plt.errorbar(x06, y06, yerr=yerr06, xerr=xerr06, fmt='^', label='2006')\n",
    "plt.errorbar(x091, y091, yerr=yerr091, xerr=xerr091, fmt='s',label='2009(1)')\n",
    "plt.errorbar(x092, y092, yerr=yerr092, xerr=xerr092, fmt='*',label='2009(2)')\n",
    "plt.xlabel(\"Energy (keV)\")\n",
    "plt.ylabel(r\"$keV^2 (Photons{\\ }cm^{-2}{\\ } s^{-1}{\\ } keV^{-1})$\")\n",
    "plt.legend()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
